(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 7.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[    122486,       3119]
NotebookOptionsPosition[    119861,       3036]
NotebookOutlinePosition[    120387,       3054]
CellTagsIndexPosition[    120344,       3051]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell["\<\
The order of variables:
{Ro, polyp, polypi, NS5B, NS4B, NS4NS5, complex, cellfactor, complexwork, \
Nfull, Nempty, NS4NS5adsfast, NS4NS5adsslow, ingNS3}\
\>", "Text",
 CellChangeTimes->{{3.5086702825*^9, 3.508670282625*^9}, {3.508670377625*^9, 
   3.508670400421875*^9}, {3.50867065784375*^9, 3.508670667578125*^9}, {
   3.50867072834375*^9, 3.508670812640625*^9}, {3.508670956125*^9, 
   3.50867108778125*^9}, {3.50867150575*^9, 3.50867150675*^9}, 
   3.525441877935867*^9}],

Cell["Version :: 'equations 11_05_2011_ ispravl versiya.docx':", "Text",
 CellChangeTimes->{{3.515490093171875*^9, 3.515490132625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"variables", " ", "in", " ", "static", " ", 
    RowBox[{"form", " ", "--"}], " ", "for", " ", "fixed", " ", "point", " ", 
    RowBox[{"calculation", ":"}]}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"var", "=", 
     RowBox[{"{", 
      RowBox[{
      "Ro", ",", "polyp", ",", "polypi", ",", "NS5B", ",", "NS4B", ",", 
       "NS4NS5", ",", "complex", ",", "cellfactor", ",", "complexwork", ",", 
       "Nfull", ",", "Nempty", ",", "NS4NS5adsfast", ",", "NS4NS5adsslow", 
       ",", "ingNS3"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "RHS", " ", "of", " ", "the", " ", "equations", " ", "in", " ", "static", 
     " ", 
     RowBox[{"form", ":"}]}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"F", "=", 
     RowBox[{"Table", "[", 
      RowBox[{"0", ",", 
       RowBox[{"{", "14", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"F", "[", 
      RowBox[{"[", "1", "]"}], "]"}], "=", 
     RowBox[{
      RowBox[{"kout", "*", "Nfull"}], "-", 
      RowBox[{"krnaads", "*", "complex", "*", "Ro"}], "+", 
      RowBox[{"krnadeads", "*", "complexwork"}], "-", 
      RowBox[{"mr", "*", "Ro"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"F", "[", 
      RowBox[{"[", "2", "]"}], "]"}], "=", 
     RowBox[{
      RowBox[{"k2", "*", "Ro"}], "+", 
      RowBox[{"kpolypiobr", "*", "polypi"}], "-", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"kc", "+", "mpolyp"}], ")"}], "*", "polyp"}], "-", 
      RowBox[{"kpolypi", "*", "ingNS3", "*", "polyp"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"F", "[", 
      RowBox[{"[", "3", "]"}], "]"}], "=", 
     RowBox[{
      RowBox[{"kpolypi", "*", "ingNS3", "*", "polyp"}], "-", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"kpolypiobr", "+", "mpolypi"}], ")"}], "*", "polypi"}]}]}], 
    ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"F", "[", 
      RowBox[{"[", "4", "]"}], "]"}], "=", 
     RowBox[{
      RowBox[{"kc", "*", "polyp"}], "-", 
      RowBox[{"vesNS5B", "*", "kcomplex", "*", "NS4B", "*", "NS5B"}], "-", 
      RowBox[{"me", "*", "NS5B"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"F", "[", 
      RowBox[{"[", "5", "]"}], "]"}], "=", 
     RowBox[{
      RowBox[{"kc", "*", "polyp"}], "-", 
      RowBox[{"vesNS4B", "*", "kcomplex", "*", "NS4B", "*", "NS5B"}], "-", 
      RowBox[{"me", "*", "NS4B"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"F", "[", 
      RowBox[{"[", "6", "]"}], "]"}], "=", 
     RowBox[{
      RowBox[{"kcomplex", "*", "NS4B", "*", "NS5B"}], "-", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"kadsfast", "+", "me"}], ")"}], "*", "NS4NS5"}], "-", 
      RowBox[{"kads", "*", "NS4NS5", "*", "cellfactor"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"F", "[", 
      RowBox[{"[", "7", "]"}], "]"}], "=", 
     RowBox[{
      RowBox[{"kads", "*", "NS4NS5", "*", "cellfactor"}], "-", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"kv", "+", "me"}], ")"}], "*", "complex"}], "+", 
      RowBox[{"kdeads", "*", "cellfactor", "*", 
       RowBox[{"(", 
        RowBox[{"NS4NS5adsfast", "+", "NS4NS5adsslow"}], ")"}]}], "-", 
      RowBox[{"krnaads", "*", "complex", "*", "Ro"}], "+", 
      RowBox[{"krnadeads", "*", "complexwork"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"F", "[", 
      RowBox[{"[", "8", "]"}], "]"}], "=", 
     RowBox[{"kcfin", "-", 
      RowBox[{"kads", "*", "NS4NS5", "*", "cellfactor"}], "-", 
      RowBox[{"kdeads", "*", "cellfactor", "*", 
       RowBox[{"(", 
        RowBox[{"NS4NS5adsfast", "+", "NS4NS5adsslow"}], ")"}]}], "-", 
      RowBox[{"mcf", "*", "cellfactor"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"F", "[", 
      RowBox[{"[", "9", "]"}], "]"}], "=", 
     RowBox[{
      RowBox[{"krnaads", "*", "complex", "*", "Ro"}], "-", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"krnadeads", "+", "kv", "+", "me"}], ")"}], "*", 
       "complexwork"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"F", "[", 
      RowBox[{"[", "10", "]"}], "]"}], "=", 
     RowBox[{
      RowBox[{"kv", "*", "complexwork"}], "-", 
      RowBox[{"me", "*", "Nfull"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"F", "[", 
      RowBox[{"[", "11", "]"}], "]"}], "=", 
     RowBox[{
      RowBox[{"kv", "*", "complex"}], "-", 
      RowBox[{"me", "*", "Nempty"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"F", "[", 
      RowBox[{"[", "12", "]"}], "]"}], "=", 
     RowBox[{
      RowBox[{"kcheck", "*", "kdiff", "*", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "NS4NS5adsfast"}], "+", "NS4NS5adsslow"}], ")"}]}], "+", 
      RowBox[{"kadsfast", "*", "NS4NS5"}], "-", 
      RowBox[{"kdeads", "*", "cellfactor", "*", "NS4NS5adsfast"}], "-", 
      RowBox[{"me", "*", "NS4NS5adsfast"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"F", "[", 
      RowBox[{"[", "13", "]"}], "]"}], "=", 
     RowBox[{"kcheck", "*", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"kdiff", "*", 
         RowBox[{"(", 
          RowBox[{"NS4NS5adsfast", "-", "NS4NS5adsslow"}], ")"}]}], "-", 
        RowBox[{"kdeads", "*", "cellfactor", "*", "NS4NS5adsslow"}], "-", 
        RowBox[{"mlipidslow", "*", "NS4NS5adsslow"}]}], ")"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"F", "[", 
      RowBox[{"[", "14", "]"}], "]"}], "=", 
     RowBox[{"kdrugin", "-", 
      RowBox[{"mdrug", "*", "ingNS3"}], "-", 
      RowBox[{"kpolypi", "*", "ingNS3", "*", "polyp"}], "+", 
      RowBox[{"kpolypiobr", "*", "polypi"}]}]}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.5086705564375*^9, 3.508670655453125*^9}, {
   3.508671112984375*^9, 3.5086712325625*^9}, {3.50867127715625*^9, 
   3.50867131075*^9}, {3.508671424421875*^9, 3.50867173209375*^9}, {
   3.508671794140625*^9, 3.508671846546875*^9}, {3.508671890765625*^9, 
   3.5086720038125*^9}, {3.508672051234375*^9, 3.50867218634375*^9}, {
   3.508672218765625*^9, 3.508672265453125*^9}, {3.5086723639375*^9, 
   3.508672406828125*^9}, {3.508672442203125*^9, 3.508672639375*^9}, {
   3.5086726735625*^9, 3.508672686078125*^9}, 3.5088291941875*^9, {
   3.508829618640625*^9, 3.5088302688125*^9}, {3.50883030034375*^9, 
   3.50883030575*^9}, 3.511426717671875*^9, {3.511430847359375*^9, 
   3.51143090009375*^9}, {3.5115155280307484`*^9, 3.511515560728406*^9}, {
   3.5115156579633765`*^9, 3.51151566556059*^9}, {3.5115220600560217`*^9, 
   3.511522069618839*^9}, {3.511545786974783*^9, 3.5115457897919445`*^9}, {
   3.511546029784671*^9, 3.5115460327348394`*^9}, 3.515490064265625*^9, {
   3.515490340375*^9, 3.51549040121875*^9}, 3.51549055753125*^9}],

Cell["\<\
Parameters fixed from literature :: k2, me, kpolypi, kpolypi_obr, kc, \
kdrug_in.\
\>", "Text",
 CellChangeTimes->{{3.515933681796875*^9, 3.515933767890625*^9}, {
   3.515933887578125*^9, 3.515933894625*^9}, 3.5160967770625*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"parameters", " ", "fixed", " ", "from", " ", 
    RowBox[{"literature", ":"}]}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"k2", ":=", "150"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"kc", ":=", "1.0"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"me", ":=", "0.08"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"kpolypi", ":=", "0.0000307"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"kpolypiobr", ":=", 
     RowBox[{"42000", "*", "kpolypi"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"kdrugin", ":=", 
     RowBox[{"mdrug", "*", "ingNS3conc"}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"parameters", " ", "already", " ", 
     RowBox[{"studied", ":"}]}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"kcomplex", ":=", "25.1"}], ";"}], 
   RowBox[{"(*", " ", 
    RowBox[{
    "leave", " ", "the", " ", "optimal", " ", "value", " ", "since", " ", 
     "there", " ", "is", " ", "no", " ", "dependence", " ", "on", " ", "it"}],
     " ", "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"parameters", " ", "dependent", " ", "on", " ", "other", " ", 
     RowBox[{"parameters", ":"}]}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"kads", ":=", 
     RowBox[{"kadscoef", "*", "kadsfast"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"kcfin", ":=", 
     RowBox[{"kcfincoef", "*", "mcf", "*", "interferoncf"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"mr", ":=", 
     RowBox[{"mrcoef", "*", "interferonmr"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"vesNS5B", ":=", "vesNS4B"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"interferonmr", ":=", 
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"kdrugin", "\[Equal]", "0"}], ",", "1", ",", 
       RowBox[{"interferonmrA", "+", 
        RowBox[{"interferonmrB", "/", 
         RowBox[{"(", 
          RowBox[{"interferonmrC", "+", 
           RowBox[{"Exp", "[", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"polyp", "-", "interferonmrD"}], ")"}], "/", 
             "interferonmrE"}], "]"}]}], ")"}]}]}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"interferoncf", ":=", 
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"kdrugin", "\[Equal]", "0"}], ",", "1", ",", 
       RowBox[{"interferoncfA", "+", 
        RowBox[{"interferoncfB", "/", 
         RowBox[{"(", 
          RowBox[{"interferoncfC", "+", 
           RowBox[{"Exp", "[", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"polyp", "-", "interferoncfD"}], ")"}], "/", 
             "interferoncfE"}], "]"}]}], ")"}]}]}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
     RowBox[{"parameters", " ", "to", " ", 
      RowBox[{"vary", " ", ":", " ", 
       RowBox[{
       "some", " ", "parameters", " ", "depend", " ", "on", " ", 
        "others"}]}]}], ",", " ", 
     RowBox[{
     "so", " ", "we", " ", "vary", " ", "only", " ", "corresponding", " ", 
      "multipliers", " ", "for", " ", "these", " ", 
      RowBox[{"parameters", ":"}]}]}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"params", "=", 
     RowBox[{"{", 
      RowBox[{
      "kout", ",", "kv", ",", "krnaads", ",", "krnadeads", ",", "kadsfast", 
       ",", "kdeads", ",", "mlipidslow", ",", "kadscoef", ",", "mcf", ",", 
       "kcfincoef", ",", "mpolyp", ",", "mpolypi", ",", "mrcoef", ",", 
       "vesNS4B", ",", "kdiff", ",", "mdrug", ",", "interferonmrA", ",", 
       "interferonmrB", ",", "interferonmrC", ",", "interferonmrD", ",", 
       "interferonmrE", ",", "interferoncfA", ",", "interferoncfB", ",", 
       "interferoncfC", ",", "interferoncfD", ",", "interferoncfE"}], "}"}]}],
     ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"optimal", " ", 
     RowBox[{"values", ":"}]}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"optval", "=", 
     RowBox[{"{", 
      RowBox[{"4", ",", "1.5", ",", "4", ",", "0.1", ",", "20.12", ",", 
       RowBox[{"5.8909", "*", 
        SuperscriptBox["10", 
         RowBox[{"-", "5"}]]}], ",", "0.001", ",", "0.1", ",", "0.4792", ",", 
       "18.7829", ",", "0", ",", "0.2778", ",", "1", ",", "1000", ",", 
       RowBox[{"1.3333", "*", 
        SuperscriptBox["10", 
         RowBox[{"-", "4"}]]}], ",", "0.1", ",", "1", ",", "0.5", ",", "1", 
       ",", "8.54", ",", "8.19", ",", "1", ",", 
       RowBox[{"-", "1"}], ",", "1", ",", "8.54", ",", "8.19"}], "}"}]}], 
    ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"kcheck", "=", "1"}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.516096821984375*^9, 3.516096837203125*^9}, {
   3.5160968906875*^9, 3.5160969810625*^9}, {3.516097135546875*^9, 
   3.51609716496875*^9}, {3.516097267609375*^9, 3.51609732478125*^9}, {
   3.516097704890625*^9, 3.516097781578125*^9}, {3.5160978159375*^9, 
   3.516097819734375*^9}, {3.516098340765625*^9, 3.51609841475*^9}, {
   3.51609846671875*^9, 3.51609850465625*^9}, {3.51609887115625*^9, 
   3.51609888428125*^9}, {3.516098940984375*^9, 3.516099098640625*^9}, {
   3.516099278015625*^9, 3.51609936621875*^9}, {3.516099761625*^9, 
   3.516099828484375*^9}, {3.5160998615*^9, 3.516099888953125*^9}, {
   3.51610001934375*^9, 3.516100021140625*^9}, {3.51610037640625*^9, 
   3.516100378515625*^9}, {3.51610238765625*^9, 3.516102481640625*^9}, {
   3.516102511671875*^9, 3.51610253109375*^9}, {3.516102568484375*^9, 
   3.516102576828125*^9}, {3.516102653703125*^9, 3.516102673265625*^9}, {
   3.51610272821875*^9, 3.51610280465625*^9}, {3.516102836984375*^9, 
   3.51610294896875*^9}, {3.516103045703125*^9, 3.516103194953125*^9}, {
   3.51610325225*^9, 3.5161034228125*^9}, {3.516105565109375*^9, 
   3.5161055694375*^9}, 3.5163552355*^9, {3.516441272140625*^9, 
   3.516441389765625*^9}, {3.516441427203125*^9, 3.516441437140625*^9}}],

Cell["RHS in dynamical form (to calculate solution dynamics):", "Text",
 CellChangeTimes->{{3.509196964854842*^9, 3.5091969917648897`*^9}, {
  3.509201151761737*^9, 3.5092011526509385`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"dynrule", ":=", 
   RowBox[{
    RowBox[{
     RowBox[{"#", "\[Rule]", 
      RowBox[{"#", "[", "t", "]"}]}], "&"}], "/@", "var"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"vardyn", ":=", 
   RowBox[{"var", "/.", "dynrule"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Fdyn", ":=", 
    RowBox[{"F", "/.", "dynrule"}]}], ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{"total", " ", "replicon", " ", "RNA", " ", 
    RowBox[{"concentration", ":"}]}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Rtotal", "[", "t_", "]"}], ":=", 
   RowBox[{
    RowBox[{"vardyn", "[", 
     RowBox[{"[", "1", "]"}], "]"}], "+", 
    RowBox[{"6", 
     RowBox[{"vardyn", "[", 
      RowBox[{"[", "10", "]"}], "]"}]}], "+", 
    RowBox[{"vardyn", "[", 
     RowBox[{"[", "9", "]"}], "]"}]}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.509197006756516*^9, 3.509197064180217*^9}, {
   3.509197097642276*^9, 3.5091971241779222`*^9}, {3.5092006138727922`*^9, 
   3.509200614371993*^9}, {3.51575959446875*^9, 3.51575959834375*^9}, {
   3.5159214939375*^9, 3.515921525109375*^9}, 3.515921791703125*^9, 
   3.516103790234375*^9}],

Cell["\<\
Calculate solution for optimal parameter values in 4 situations (no drug and \
3 drug concentrations):\
\>", "Text",
 CellChangeTimes->{{3.51610633040625*^9, 3.516106391796875*^9}, {
  3.524211269953125*^9, 3.52421129259375*^9}, {3.524212879171875*^9, 
  3.524212893125*^9}, {3.5254419555772033`*^9, 3.525441964828019*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
   "solve", " ", "for", " ", "optimal", " ", "parameter", " ", "values", " ", 
    "for", " ", "further", " ", 
    RowBox[{"comparison", ":"}]}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"tmax", "=", 
     RowBox[{"{", 
      RowBox[{"10000", ",", "365", ",", "365", ",", "365"}], "}"}]}], ";"}], 
   " ", 
   RowBox[{"(*", " ", 
    RowBox[{
    "the", " ", "drug", " ", "treatment", " ", "time", " ", "is", " ", "362", 
     RowBox[{"h", "."}]}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Do", "[", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{
       RowBox[{"ingNS3conc", "=", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{"0", ",", 
           RowBox[{"1.5", "*", 
            SuperscriptBox["10", "6"]}], ",", 
           RowBox[{"7.5", "*", 
            SuperscriptBox["10", "6"]}], ",", 
           RowBox[{"15", "*", 
            SuperscriptBox["10", "6"]}]}], "}"}], "[", 
         RowBox[{"[", "n", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"Fdyn2", "=", 
        RowBox[{"Fdyn", "/.", 
         RowBox[{"Table", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"params", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "\[Rule]", 
            RowBox[{"optval", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"i", ",", 
             RowBox[{"Length", "@", "params"}]}], "}"}]}], "]"}]}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"incond", "=", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{"n", "\[Equal]", "1"}], ",", "\[IndentingNewLine]", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{
             RowBox[{"Ro", "[", "0", "]"}], "\[Equal]", "10"}], ",", 
            RowBox[{
             RowBox[{"polyp", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
            RowBox[{
             RowBox[{"polypi", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
            RowBox[{
             RowBox[{"NS5B", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
            RowBox[{
             RowBox[{"NS4B", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
            RowBox[{
             RowBox[{"NS4NS5", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
            RowBox[{
             RowBox[{"complex", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
            RowBox[{
             RowBox[{"cellfactor", "[", "0", "]"}], "\[Equal]", "18"}], ",", 
            RowBox[{
             RowBox[{"complexwork", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
            RowBox[{
             RowBox[{"Nfull", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
            RowBox[{
             RowBox[{"Nempty", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
            RowBox[{
             RowBox[{"NS4NS5adsfast", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
            RowBox[{
             RowBox[{"NS4NS5adsslow", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
            RowBox[{
             RowBox[{"ingNS3", "[", "0", "]"}], "\[Equal]", "0"}]}], "}"}], 
          ",", "\[IndentingNewLine]", 
          RowBox[{"Table", "[", 
           RowBox[{
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{"vardyn", "/.", 
                RowBox[{"t", "\[Rule]", "0"}]}], ")"}], "[", 
              RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"vardyn", "/.", 
                 RowBox[{"sol", "[", "1", "]"}]}], "/.", 
                RowBox[{"t", "\[Rule]", 
                 RowBox[{"tmax", "[", 
                  RowBox[{"[", "1", "]"}], "]"}]}]}], ")"}], "[", 
              RowBox[{"[", "i", "]"}], "]"}]}], ",", 
            RowBox[{"{", 
             RowBox[{"i", ",", "14"}], "}"}]}], "]"}]}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"eqs", "=", 
        RowBox[{
         RowBox[{"Table", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"D", "[", 
             RowBox[{
              RowBox[{"vardyn", "[", 
               RowBox[{"[", "i", "]"}], "]"}], ",", "t"}], "]"}], "\[Equal]", 
            RowBox[{"Fdyn2", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"i", ",", 
             RowBox[{"Length", "@", "vardyn"}]}], "}"}]}], "]"}], "~", "Join",
          "~", "incond"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"sol", "[", "n", "]"}], "=", 
        RowBox[{"First", "@", 
         RowBox[{"NDSolve", "[", 
          RowBox[{"eqs", ",", "vardyn", ",", 
           RowBox[{"{", 
            RowBox[{"t", ",", "0", ",", 
             RowBox[{"tmax", "[", 
              RowBox[{"[", "n", "]"}], "]"}]}], "}"}], ",", 
           RowBox[{"InterpolationOrder", "\[Rule]", "1"}]}], "]"}]}]}]}], ",",
       "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{"n", ",", "4"}], "}"}]}], "]"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
     RowBox[{"steady", " ", "state", " ", "for", " ", "no"}], "-", 
     RowBox[{"drug", " ", 
      RowBox[{"solution", ":"}]}]}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"RtCheck1", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"Rtotal", "[", "t", "]"}], "/.", 
       RowBox[{"sol", "[", "1", "]"}]}], "/.", 
      RowBox[{"t", "\[Rule]", 
       RowBox[{"tmax", "[", 
        RowBox[{"[", "1", "]"}], "]"}]}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "and", " ", "53", " ", "control", " ", "time", " ", "points", " ", "for", 
     " ", "solutions", " ", "under", " ", "drugs", " ", "from", " ", "the", 
     " ", "data", " ", "time", " ", 
     RowBox[{"period", ":"}]}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"RtCheck2", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"Rtotal", "[", "t", "]"}], "/.", 
         RowBox[{"sol", "[", "n", "]"}]}], "/.", 
        RowBox[{"t", "\[Rule]", "tt"}]}], ",", 
       RowBox[{"{", 
        RowBox[{"n", ",", "2", ",", "4", ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"tt", ",", "0", ",", "365", ",", "7"}], "}"}]}], "]"}]}], 
    ";"}]}]}]], "Input",
 CellChangeTimes->{{3.51154052893004*^9, 3.5115405454769864`*^9}, {
   3.5115433063859015`*^9, 3.511543342163948*^9}, {3.5115465004895935`*^9, 
   3.511546503992794*^9}, {3.5115465882206116`*^9, 3.511546590759757*^9}, {
   3.511598239140625*^9, 3.511598242203125*^9}, {3.515491826609375*^9, 
   3.515491826875*^9}, {3.516103839171875*^9, 3.51610395640625*^9}, {
   3.516104024109375*^9, 3.516104048375*^9}, {3.516104091921875*^9, 
   3.516104307046875*^9}, {3.51610455890625*^9, 3.51610459965625*^9}, {
   3.5161047974375*^9, 3.516104820765625*^9}, {3.51610507159375*^9, 
   3.51610511646875*^9}, 3.516105195890625*^9, 3.516105426140625*^9, 
   3.51610563915625*^9, 3.5161056893125*^9, 3.51610585353125*^9, {
   3.51610606271875*^9, 3.516106091828125*^9}, {3.516353670609375*^9, 
   3.516353824140625*^9}, {3.51635386859375*^9, 3.516353896359375*^9}, {
   3.51635406440625*^9, 3.51635414846875*^9}, {3.516354268296875*^9, 
   3.5163542969375*^9}, {3.51635433075*^9, 3.5163544541875*^9}, {
   3.5166157804375*^9, 3.51661581375*^9}, {3.516615940609375*^9, 
   3.516615959*^9}, {3.524212903984375*^9, 3.524212904734375*^9}, {
   3.5242131545*^9, 3.52421315775*^9}, {3.5242141874375*^9, 
   3.524214188484375*^9}, {3.5242144653125*^9, 3.52421446615625*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{"RtCheck1", "\[IndentingNewLine]", "RtCheck2"}], "Input",
 CellChangeTimes->{{3.516355275234375*^9, 3.5163552854375*^9}}],

Cell[BoxData["969.1974452960895`"], "Output",
 CellChangeTimes->{3.51635528609375*^9, 3.516441498453125*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
    "969.1974452960895`", ",", "869.1193953940244`", ",", 
     "739.9482405147216`", ",", "618.9674102318944`", ",", 
     "512.6154525952649`", ",", "421.9578382387168`", ",", 
     "345.9661527661351`", ",", "282.94622160635544`", ",", 
     "231.0894856779158`", ",", "188.67941319078918`", ",", 
     "154.16716020544533`", ",", "126.19497256622748`", ",", 
     "103.5961807534396`", ",", "85.38359342661636`", ",", 
     "70.73212996309263`", ",", "58.959098300210364`", ",", 
     "49.50424790402063`", ",", "41.910892491057666`", ",", 
     "35.808808815361786`", ",", "30.89922342276062`", ",", 
     "26.94192151215926`", ",", "23.74436691801823`", ",", 
     "21.152630916602206`", ",", "19.04390464307136`", ",", 
     "17.32036498435412`", ",", "15.904177765333404`", ",", 
     "14.733451164348786`", ",", "13.758976764029935`", ",", 
     "12.941618383664732`", ",", "12.250231813066167`", ",", 
     "11.660011512218558`", ",", "11.151183604495195`", ",", 
     "10.707969189094953`", ",", "10.317764248095113`", ",", 
     "9.970489663492286`", ",", "9.658075941008939`", ",", 
     "9.374052704085743`", ",", "9.113224771347571`", ",", 
     "8.871412403706294`", ",", "8.645246044409358`", ",", 
     "8.432002643714188`", ",", "8.229474781017217`", ",", 
     "8.035866655095836`", ",", "7.84971109705687`", ",", 
     "7.6698028297925145`", ",", "7.495145297049898`", ",", 
     "7.32490784756069`", ",", "7.158391296599979`", ",", 
     "6.9949999140196395`", ",", "6.834218850196461`", ",", 
     "6.675595544383148`", ",", "6.518724338267809`", ",", 
     "6.363233628550176`", ",", "6.208774881100206`", ",", 
     "6.055013039620971`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "969.1974452960895`", ",", "715.9481186980696`", ",", 
     "477.77743201038754`", ",", "314.2818850463402`", ",", 
     "205.81379207485185`", ",", "135.06375663062266`", ",", 
     "89.3430581989227`", ",", "59.9537448999578`", ",", "41.11256601496832`",
      ",", "29.01474761001299`", ",", "21.13312168810559`", ",", 
     "15.807029033265916`", ",", "12.023289386139165`", ",", 
     "9.218674785916527`", ",", "7.093194126906077`", ",", 
     "5.478423511965871`", ",", "4.261261237986458`", ",", 
     "3.3514489564939103`", ",", "2.6739653513591337`", ",", 
     "2.168476336615522`", ",", "1.7885669949320637`", ",", 
     "1.499781124303214`", ",", "1.277141473108997`", ",", 
     "1.1027794586000739`", ",", "0.9639732829225822`", ",", 
     "0.8516524671696002`", ",", "0.7593121892178893`", ",", 
     "0.6822469726544167`", ",", "0.6170173491632088`", ",", 
     "0.5610808001689469`", ",", "0.5125359936709575`", ",", 
     "0.46994475303344097`", ",", "0.4322070939051619`", ",", 
     "0.39847282781336857`", ",", "0.3680782516149008`", ",", 
     "0.34050037963547286`", ",", "0.31532355400116174`", ",", 
     "0.2922147357743229`", ",", "0.27090506941351244`", ",", 
     "0.2511759909954029`", ",", "0.2328486363227063`", ",", 
     "0.21577570618938982`", ",", "0.19983514184704004`", ",", 
     "0.1849251810279663`", ",", "0.17096045248243597`", ",", 
     "0.15786885516503035`", ",", "0.14558906238861885`", ",", 
     "0.13406849764132672`", ",", "0.12326168417077402`", ",", 
     "0.11312889052295744`", ",", "0.10363501615765104`", ",", 
     "0.09474865875962662`", ",", "0.08644135807085743`", ",", 
     "0.07868694589029654`", ",", "0.07146102308666956`"}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.51635528609375*^9, 3.51644149846875*^9}]
}, Open  ]],

Cell["Prescribed ranges:", "Text",
 CellChangeTimes->{{3.5242158156875*^9, 3.524215824953125*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"ranges", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"3.8034806490", ",", "4.1965193510"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"1.1500978957", ",", "1.8499021043"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"3.4201974172", ",", "4.5798025828"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.0805358509", ",", "0.1194641491"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"12.6454302489", ",", "27.5945697511"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.0000545067", ",", "0.0000633133"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.0005000000", ",", "0.0233410738"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.0908145032", ",", "0.1091854968"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.4666796563", ",", "0.4917203437"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"18.3357287261", ",", "19.2300712739"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.0000000000", ",", "0.0674968579"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.2761456965", ",", "0.2794543035"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.9909333282", ",", "1.0090666718"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"991.0617287694", ",", "1008.9382712306"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.0000666500", ",", "0.0005881878"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.0926216595", ",", "0.1073783405"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.9951943629", ",", "1.0048056371"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.4494649602", ",", "0.5505350398"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.5188481930", ",", "1.4811518070"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"8.2385836414", ",", "8.8414163586"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"7.1340472810", ",", "9.2459527190"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.9050401809", ",", "1.0949598191"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "1.2319611993"}], ",", 
       RowBox[{"-", "0.7680388007"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0.8457806932", ",", "1.1542193068"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"4.2700000000", ",", "14.3480146623"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"7.4087580297", ",", "8.9712419703"}], "}"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"rangesInit", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"9.4453688487", ",", "10.5546311513"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"12.3890986149", ",", "23.6109013851"}], "}"}]}], "}"}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.524216242890625*^9, 3.52421629915625*^9}, {
  3.524216336578125*^9, 3.52421658053125*^9}, {3.52421675053125*^9, 
  3.52421698346875*^9}, {3.52421701978125*^9, 3.524217154328125*^9}, {
  3.5248930948125*^9, 3.524893354453125*^9}, {3.52489358821875*^9, 
  3.524893707046875*^9}, {3.52489377525*^9, 3.524893821203125*^9}, {
  3.52489388321875*^9, 3.524893900140625*^9}, {3.524893972921875*^9, 
  3.52489398284375*^9}, {3.524894026140625*^9, 3.5248940285*^9}, {
  3.5248940926875*^9, 3.524894120421875*^9}, {3.524894671890625*^9, 
  3.524894682375*^9}, {3.5248947650625*^9, 3.524894798875*^9}}],

Cell["1000 random parameter values from the ranges:", "Text",
 CellChangeTimes->{
  3.5167024455625*^9, {3.516702572734375*^9, 3.516702673484375*^9}, {
   3.516703048984375*^9, 3.516703072609375*^9}, {3.516703539703125*^9, 
   3.516703540703125*^9}, {3.51670368184375*^9, 3.516703693078125*^9}, {
   3.5168784194907007`*^9, 3.516878423831949*^9}, {3.516878466735403*^9, 
   3.5168785349663057`*^9}, {3.5168786170129986`*^9, 3.516878636658122*^9}, {
   3.524217579796875*^9, 3.524217638203125*^9}, {3.5242925805*^9, 
   3.524292581109375*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"SeedRandom", "[", "56", "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Num", "=", "1000"}], ";", 
  RowBox[{"(*", " ", 
   RowBox[{"number", " ", "of", " ", "samples"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"tmax", "=", 
   RowBox[{"{", 
    RowBox[{"10000", ",", "365", ",", "365", ",", "365"}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"icBasic", ":=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"Ro", "[", "0", "]"}], "\[Equal]", "RoInit"}], ",", 
     RowBox[{
      RowBox[{"polyp", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
     RowBox[{
      RowBox[{"polypi", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
     RowBox[{
      RowBox[{"NS5B", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
     RowBox[{
      RowBox[{"NS4B", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
     RowBox[{
      RowBox[{"NS4NS5", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
     RowBox[{
      RowBox[{"complex", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
     RowBox[{
      RowBox[{"cellfactor", "[", "0", "]"}], "\[Equal]", "cellInit"}], ",", 
     RowBox[{
      RowBox[{"complexwork", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
     RowBox[{
      RowBox[{"Nfull", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
     RowBox[{
      RowBox[{"Nempty", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
     RowBox[{
      RowBox[{"NS4NS5adsfast", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
     RowBox[{
      RowBox[{"NS4NS5adsslow", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
     RowBox[{
      RowBox[{"ingNS3", "[", "0", "]"}], "\[Equal]", "0"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"count", "=", "0"}], ";", 
  RowBox[{"(*", " ", 
   RowBox[{"number", " ", "of", " ", "goodies"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"Do", "[", "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"select", " ", "random", " ", "parameter", " ", 
     RowBox[{"values", ":"}]}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"RoInit", "=", 
      RowBox[{"RandomReal", "[", 
       RowBox[{"rangesInit", "[", 
        RowBox[{"[", "1", "]"}], "]"}], "]"}]}], ";", 
     RowBox[{"cellInit", "=", 
      RowBox[{"RandomReal", "[", 
       RowBox[{"rangesInit", "[", 
        RowBox[{"[", "2", "]"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"randval", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"RandomReal", "[", 
         RowBox[{"ranges", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", 
          RowBox[{"Length", "@", "params"}]}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"(*", " ", 
      RowBox[{"write", " ", "these", " ", "values", " ", "to", " ", 
       RowBox[{"file", ":"}]}], " ", "*)"}], "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"RoInit", ",", "cellInit", ",", "randval"}], "}"}], ">>>", 
      RowBox[{
      "disk", "<>", 
       "\"\<Vitalik\\\\Hepatitis\\\\Math\\\\abs_log_diff\\\\dt\\\\par\\\\par_\
\>\"", "<>", 
       RowBox[{"ToString", "[", "j", "]"}]}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"Do", "[", 
      RowBox[{"(*", " ", 
       RowBox[{"cycle", " ", "over", " ", "drug", " ", "concentrations"}], 
       " ", "*)"}], "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{
        RowBox[{"ingNS3conc", "=", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"0", ",", 
            RowBox[{"1.5", "*", 
             SuperscriptBox["10", "6"]}], ",", 
            RowBox[{"7.5", "*", 
             SuperscriptBox["10", "6"]}], ",", 
            RowBox[{"15", "*", 
             SuperscriptBox["10", "6"]}]}], "}"}], "[", 
          RowBox[{"[", "n", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
        RowBox[{"Fdyn2", "=", 
         RowBox[{"Fdyn", "/.", 
          RowBox[{"Table", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"params", "[", 
              RowBox[{"[", "i", "]"}], "]"}], "\[Rule]", 
             RowBox[{"randval", "[", 
              RowBox[{"[", "i", "]"}], "]"}]}], ",", 
            RowBox[{"{", 
             RowBox[{"i", ",", 
              RowBox[{"Length", "@", "params"}]}], "}"}]}], "]"}]}]}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"incond", "=", 
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{"n", "\[Equal]", "1"}], ",", "\[IndentingNewLine]", 
           "icBasic", ",", "\[IndentingNewLine]", 
           RowBox[{"Table", "[", 
            RowBox[{
             RowBox[{
              RowBox[{
               RowBox[{"(", 
                RowBox[{"vardyn", "/.", 
                 RowBox[{"t", "\[Rule]", "0"}]}], ")"}], "[", 
               RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"vardyn", "/.", 
                  RowBox[{"sol", "[", "1", "]"}]}], "/.", 
                 RowBox[{"t", "\[Rule]", 
                  RowBox[{"tmax", "[", 
                   RowBox[{"[", "1", "]"}], "]"}]}]}], ")"}], "[", 
               RowBox[{"[", "i", "]"}], "]"}]}], ",", 
             RowBox[{"{", 
              RowBox[{"i", ",", "14"}], "}"}]}], "]"}]}], "]"}]}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"eqs", "=", 
         RowBox[{
          RowBox[{"Table", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"D", "[", 
              RowBox[{
               RowBox[{"vardyn", "[", 
                RowBox[{"[", "i", "]"}], "]"}], ",", "t"}], "]"}], "\[Equal]", 
             RowBox[{"Fdyn2", "[", 
              RowBox[{"[", "i", "]"}], "]"}]}], ",", 
            RowBox[{"{", 
             RowBox[{"i", ",", 
              RowBox[{"Length", "@", "vardyn"}]}], "}"}]}], "]"}], "~", 
          "Join", "~", "incond"}]}], ";", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"sol", "[", "n", "]"}], "=", 
         RowBox[{"First", "@", 
          RowBox[{"NDSolve", "[", 
           RowBox[{"eqs", ",", "vardyn", ",", 
            RowBox[{"{", 
             RowBox[{"t", ",", "0", ",", 
              RowBox[{"tmax", "[", 
               RowBox[{"[", "n", "]"}], "]"}]}], "}"}], ",", 
            RowBox[{"InterpolationOrder", "\[Rule]", "1"}]}], "]"}]}]}]}], 
       ",", "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{"n", ",", "4"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"(*", " ", 
      RowBox[{"collect", " ", "statistics", " ", "of", " ", "solution", " ", 
       RowBox[{"properties", ":"}]}], " ", "*)"}], "\[IndentingNewLine]", 
     RowBox[{"(*", " ", 
      RowBox[{"stability", " ", "without", " ", 
       RowBox[{"drugs", ":"}]}], " ", "*)"}], "\[IndentingNewLine]", 
     RowBox[{"stab", "=", 
      RowBox[{"Max", "[", 
       RowBox[{"Abs", "[", 
        RowBox[{
         RowBox[{"D", "[", 
          RowBox[{
           RowBox[{"vardyn", "/.", 
            RowBox[{"sol", "[", "1", "]"}]}], ",", "t"}], "]"}], "/.", 
         RowBox[{"t", "\[Rule]", 
          RowBox[{"tmax", "[", 
           RowBox[{"[", "1", "]"}], "]"}]}]}], "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"(*", " ", 
      RowBox[{"deviation", " ", "from", " ", "the", " ", "optimal", " ", 
       RowBox[{"solution", ":"}]}], " ", "*)"}], "\[IndentingNewLine]", 
     RowBox[{"solstat", "=", 
      RowBox[{
       RowBox[{
        RowBox[{"Rtotal", "[", "t", "]"}], "/.", 
        RowBox[{"sol", "[", "1", "]"}]}], "/.", 
       RowBox[{"t", "\[Rule]", 
        RowBox[{"tmax", "[", 
         RowBox[{"[", "1", "]"}], "]"}]}]}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"nodrug", "=", 
      RowBox[{"RtCheck1", "-", "solstat"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"(*", " ", 
      RowBox[{
       RowBox[{"compare", " ", "logs", " ", "of", " ", "ratios"}], ",", " ", 
       RowBox[{"not", " ", "the", " ", "solutions", " ", 
        RowBox[{"themselves", ":"}]}]}], " ", "*)"}], "\[IndentingNewLine]", 
     RowBox[{"drug", "=", 
      RowBox[{
       RowBox[{"Log", "[", 
        RowBox[{"10", ",", 
         FractionBox["RtCheck2", "RtCheck1"]}], "]"}], "-", 
       RowBox[{"Log", "[", 
        RowBox[{"10", ",", 
         RowBox[{"Table", "[", 
          RowBox[{
           FractionBox[
            RowBox[{
             RowBox[{
              RowBox[{"Rtotal", "[", "t", "]"}], "/.", 
              RowBox[{"sol", "[", "n", "]"}]}], "/.", 
             RowBox[{"t", "\[Rule]", "tt"}]}], "solstat"], ",", 
           RowBox[{"{", 
            RowBox[{"n", ",", "2", ",", "4", ",", "1"}], "}"}], ",", 
           RowBox[{"{", 
            RowBox[{"tt", ",", "0", ",", "365", ",", "7"}], "}"}]}], "]"}]}], 
        "]"}]}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"(*", " ", 
      RowBox[{
      "right", " ", "statistics", " ", "and", " ", "solution", " ", "at", " ",
        "tmax", " ", "to", " ", 
       RowBox[{"file", ":"}]}], " ", "*)"}], "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"stab", ",", "nodrug", ",", "drug", ",", 
        RowBox[{"Table", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"vardyn", "/.", 
            RowBox[{"sol", "[", "n", "]"}]}], "/.", 
           RowBox[{"t", "\[Rule]", 
            RowBox[{"tmax", "[", 
             RowBox[{"[", "n", "]"}], "]"}]}]}], ",", 
          RowBox[{"{", 
           RowBox[{"n", ",", "4"}], "}"}]}], "]"}]}], "}"}], ">>>", 
      RowBox[{
      "disk", "<>", 
       "\"\<Vitalik\\\\Hepatitis\\\\Math\\\\abs_log_diff\\\\dt\\\\sol\\\\sol_\
\>\"", "<>", 
       RowBox[{"ToString", "[", "j", "]"}]}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"(*", " ", 
      RowBox[{"check", " ", "if", " ", 
       RowBox[{"it", "'"}], "s", " ", "a", " ", 
       RowBox[{"goodie", ":"}]}], " ", "*)"}], "\[IndentingNewLine]", 
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"Abs", "[", "nodrug", "]"}], "\[LessEqual]", 
         RowBox[{".3", "RtCheck1"}]}], "&&", 
        RowBox[{
         RowBox[{"Max", "[", 
          RowBox[{"Abs", "[", 
           RowBox[{"Flatten", "[", "drug", "]"}], "]"}], "]"}], 
         "\[LessEqual]", ".27"}]}], ",", 
       RowBox[{"count", "++"}]}], "]"}]}], ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{"j", ",", "Num"}], "}"}]}], "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.51670381225*^9, 3.516703820484375*^9}, 
   3.5167038566875*^9, {3.516704126375*^9, 3.516704171390625*^9}, 
   3.5167042494375*^9, {3.516704477859375*^9, 3.516704574375*^9}, {
   3.51670460853125*^9, 3.51670461425*^9}, {3.516704653046875*^9, 
   3.516704655125*^9}, {3.516704745890625*^9, 3.516704812828125*^9}, 
   3.516704849796875*^9, 3.516705020015625*^9, {3.516705149375*^9, 
   3.5167051625625*^9}, {3.516705307453125*^9, 3.5167054144375*^9}, {
   3.5167055903125*^9, 3.516705708984375*^9}, {3.51670574325*^9, 
   3.516705746109375*^9}, {3.51670583053125*^9, 3.516705865703125*^9}, {
   3.516705914703125*^9, 3.516705963625*^9}, {3.51670607303125*^9, 
   3.51670608284375*^9}, {3.51670612765625*^9, 3.516706188421875*^9}, {
   3.51670622859375*^9, 3.51670635634375*^9}, {3.5167064151875*^9, 
   3.5167064841875*^9}, {3.516706538265625*^9, 3.516706668890625*^9}, {
   3.516706779109375*^9, 3.51670680121875*^9}, {3.5168786619935713`*^9, 
   3.5168786979976306`*^9}, {3.516878761658272*^9, 3.5168787750270367`*^9}, {
   3.516878909619735*^9, 3.5168789473918953`*^9}, {3.5168789989788456`*^9, 
   3.5168789998648963`*^9}, {3.5168791398609037`*^9, 3.516879213539118*^9}, {
   3.516879277404771*^9, 3.5168792983329678`*^9}, {3.5168793571713333`*^9, 
   3.516879384359888*^9}, {3.5168795625010777`*^9, 3.516879563524136*^9}, {
   3.516879599670203*^9, 3.516879732602807*^9}, {3.516880184468652*^9, 
   3.516880185905734*^9}, {3.516880437149104*^9, 3.516880510691311*^9}, {
   3.5168811555991974`*^9, 3.516881216645689*^9}, {3.524217644921875*^9, 
   3.524217740390625*^9}, {3.524217839125*^9, 3.524217933609375*^9}, {
   3.524218508546875*^9, 3.5242185879375*^9}, {3.52421902846875*^9, 
   3.52421903740625*^9}, {3.524219069875*^9, 3.524219148609375*^9}, {
   3.5242192649375*^9, 3.524219265828125*^9}, {3.524219296453125*^9, 
   3.52421931890625*^9}, {3.52421935015625*^9, 3.524219365984375*^9}, {
   3.52421940990625*^9, 3.524219423078125*^9}, {3.524219497765625*^9, 
   3.5242195133125*^9}, {3.524389860015625*^9, 3.524389913546875*^9}, {
   3.524389970578125*^9, 3.524390002140625*^9}, {3.5243901575*^9, 
   3.524390174890625*^9}, {3.52439021728125*^9, 3.524390219671875*^9}, {
   3.524894828046875*^9, 3.52489482815625*^9}, {3.525442003439087*^9, 
   3.5254420343115416`*^9}}],

Cell[BoxData["\<\"total runs = 1000\"\>"], "Print",
 CellChangeTimes->{3.5249013525625*^9}]
}, Open  ]],

Cell["Number of good results:", "Text",
 CellChangeTimes->{{3.52429272703125*^9, 3.524292733265625*^9}}],

Cell[CellGroupData[{

Cell[BoxData["count"], "Input",
 CellChangeTimes->{{3.524292688171875*^9, 3.524292688828125*^9}}],

Cell[BoxData["2"], "Output",
 CellChangeTimes->{
  3.524292689671875*^9, 3.5243979761875*^9, {3.5249018885*^9, 
   3.52490189215625*^9}}]
}, Open  ]],

Cell["Read results:", "Text",
 CellChangeTimes->{{3.51670852925*^9, 3.516708533875*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Num", "=", "1000"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"randsol", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", "}"}], ",", 
     RowBox[{"{", "Num", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Do", "[", 
   RowBox[{"(*", " ", 
    RowBox[{"read", " ", "only", " ", 
     RowBox[{"statistics", ":"}]}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"str", "=", 
      RowBox[{"OpenRead", "[", 
       RowBox[{
       "disk", "<>", 
        "\"\<Vitalik\\\\Hepatitis\\\\Math\\\\abs_log_diff\\\\dt\\\\sol\\\\sol_\
\>\"", "<>", 
        RowBox[{"ToString", "[", "j", "]"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"randsol", "[", 
       RowBox[{"[", "j", "]"}], "]"}], "=", 
      RowBox[{"Read", "[", "str", "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"Close", "[", "str", "]"}]}], ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{"j", ",", "Num"}], "}"}]}], "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.516453491203125*^9, 3.5164535023125*^9}, {
   3.516453537796875*^9, 3.516453609875*^9}, {3.51645365803125*^9, 
   3.516453664453125*^9}, {3.51645370021875*^9, 3.516453704453125*^9}, {
   3.516453838015625*^9, 3.516453908765625*^9}, {3.51645395353125*^9, 
   3.51645395421875*^9}, {3.516454009578125*^9, 3.516454010140625*^9}, {
   3.51645406321875*^9, 3.5164540639375*^9}, 3.516454117421875*^9, 
   3.516454326625*^9, 3.516454408609375*^9, {3.516454493671875*^9, 
   3.516454507390625*^9}, 3.516454663515625*^9, {3.5165881009130397`*^9, 
   3.5165881043606453`*^9}, {3.51662390990625*^9, 3.51662391715625*^9}, {
   3.51670015446875*^9, 3.516700210109375*^9}, {3.51670868715625*^9, 
   3.516708690953125*^9}, {3.516709053046875*^9, 3.516709158140625*^9}, {
   3.516709285015625*^9, 3.516709287703125*^9}, {3.51694802809375*^9, 
   3.516948029234375*^9}, {3.516948072609375*^9, 3.5169481021875*^9}, 
   3.51694823765625*^9, 3.516948315671875*^9, 3.516948526875*^9, 
   3.516948657234375*^9, {3.516949482125*^9, 3.516949511609375*^9}, {
   3.524292631640625*^9, 3.524292666234375*^9}, {3.5243836091875*^9, 
   3.524383613640625*^9}}],

Cell["Stability:", "Text",
 CellChangeTimes->{{3.516709323375*^9, 3.51670932575*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Max", "@", 
  RowBox[{"randsol", "[", 
   RowBox[{"[", 
    RowBox[{"All", ",", "1"}], "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.516709329515625*^9, 3.516709343875*^9}}],

Cell[BoxData["0.00003383660411708083`"], "Output",
 CellChangeTimes->{3.516713359328125*^9, 3.516876073955544*^9, 
  3.51695026575*^9, 3.524304580265625*^9, 3.524398015421875*^9, 
  3.524902035828125*^9}]
}, Open  ]],

Cell["How many fulfil criteria:", "Text",
 CellChangeTimes->{{3.5167092978125*^9, 3.516709308609375*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Select", "[", 
   RowBox[{"randsol", ",", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"Abs", "[", 
        RowBox[{"#", "[", 
         RowBox[{"[", "2", "]"}], "]"}], "]"}], "\[LessEqual]", 
       RowBox[{".3", "RtCheck1"}]}], "&&", 
      RowBox[{
       RowBox[{"Max", "[", 
        RowBox[{"Abs", "[", 
         RowBox[{"Flatten", "[", 
          RowBox[{"#", "[", 
           RowBox[{"[", "3", "]"}], "]"}], "]"}], "]"}], "]"}], 
       "\[LessEqual]", ".27"}]}], "&"}]}], "]"}], "//", "Length"}]], "Input",
 CellChangeTimes->{{3.516709385875*^9, 3.51670938765625*^9}, {
  3.516710040734375*^9, 3.516710073640625*^9}, {3.51671364178125*^9, 
  3.516713643984375*^9}, {3.516714149203125*^9, 3.516714154515625*^9}, {
  3.516714190265625*^9, 3.516714209109375*^9}, {3.516876335138483*^9, 
  3.5168763826391993`*^9}, {3.5168764402084923`*^9, 3.5168764427106357`*^9}, {
  3.5168765104465103`*^9, 3.516876512173609*^9}}],

Cell[BoxData["2"], "Output",
 CellChangeTimes->{3.5168764436856914`*^9, 3.5168765130936613`*^9, 
  3.516950281078125*^9, 3.524304610953125*^9, 3.52439802121875*^9, 
  3.524902038796875*^9}]
}, Open  ]],

Cell["Distribution of no-drug steady RNA:", "Text",
 CellChangeTimes->{{3.516876648613413*^9, 3.516876660539095*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"fig1", "=", 
  RowBox[{"Histogram", "[", 
   RowBox[{
    RowBox[{"RtCheck1", "-", 
     RowBox[{"randsol", "[", 
      RowBox[{"[", 
       RowBox[{"All", ",", "2"}], "]"}], "]"}]}], ",", 
    RowBox[{"FrameLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"\"\<Steady RNA\>\"", ",", "\"\<Frequency\>\""}], "}"}]}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.516876664801338*^9, 3.5168766905718126`*^9}, {
  3.51719620278125*^9, 3.517196229359375*^9}, {3.517196508125*^9, 
  3.51719651084375*^9}}],

Cell[BoxData[
 GraphicsBox[
  {RGBColor[0.798413061722744, 0.824719615472648, 0.968322270542458], 
   EdgeForm[Opacity[0.7]], {}, 
   {RGBColor[0.798413061722744, 0.824719615472648, 0.968322270542458], 
    EdgeForm[Opacity[0.7]], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{880., 0.}, {890., 1.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{64.71312499999976, 96.63374999999996}, {
         114.90242314426742`, 118.38926705039233`}}],
       StatusArea[#, 1.]& ,
       TagBoxNote->"1."],
      StyleBox["1.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[1., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{890., 0.}, {900., 7.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{95.07124999999996, 126.99187499999971`}, {
         103.35635970751795`, 118.38926705039233`}}],
       StatusArea[#, 7.]& ,
       TagBoxNote->"7."],
      StyleBox["7.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[7., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{900., 0.}, {910., 16.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{125.42937499999971`, 157.3499999999999}, {
         86.03726455239374, 118.38926705039233`}}],
       StatusArea[#, 16.]& ,
       TagBoxNote->"16."],
      StyleBox["16.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[16., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{910., 0.}, {920., 30.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{155.7874999999999, 187.70812499999965`}, {
         59.09644986664496, 118.38926705039233`}}],
       StatusArea[#, 30.]& ,
       TagBoxNote->"30."],
      StyleBox["30.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[30., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{920., 0.}, {930., 54.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{186.14562499999965`, 218.06624999999985`}, {
         12.91219611964705, 118.38926705039233`}}],
       StatusArea[#, 54.]& ,
       TagBoxNote->"54."],
      StyleBox["54.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[54., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{930., 0.}, {940., 95.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{216.50374999999985`, 
         248.42437500000005`}, {-65.98590403147435, 118.38926705039233`}}],
       StatusArea[#, 95.]& ,
       TagBoxNote->"95."],
      StyleBox["95.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[95., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{940., 0.}, {950., 90.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{246.86187500000005`, 
         278.7824999999998}, {-56.3641845008498, 118.38926705039233`}}],
       StatusArea[#, 90.]& ,
       TagBoxNote->"90."],
      StyleBox["90.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[90., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{950., 0.}, {960., 143.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{277.2199999999998, 
         309.140625}, {-158.35441152547014`, 118.38926705039233`}}],
       StatusArea[#, 143.]& ,
       TagBoxNote->"143."],
      StyleBox["143.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[143., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{960., 0.}, {970., 138.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{307.578125, 
         339.49874999999975`}, {-148.7326919948456, 118.38926705039233`}}],
       StatusArea[#, 138.]& ,
       TagBoxNote->"138."],
      StyleBox["138.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[138., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{970., 0.}, {980., 138.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{337.93624999999975`, 
         369.85687499999995`}, {-148.7326919948456, 118.38926705039233`}}],
       StatusArea[#, 138.]& ,
       TagBoxNote->"138."],
      StyleBox["138.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[138., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{980., 0.}, {990., 120.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{368.29437499999995`, 
         400.2149999999997}, {-114.09450168459716`, 118.38926705039233`}}],
       StatusArea[#, 120.]& ,
       TagBoxNote->"120."],
      StyleBox["120.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[120., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{990., 0.}, {1000., 79.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{398.6524999999997, 
         430.5731249999999}, {-35.19640153347575, 118.38926705039233`}}],
       StatusArea[#, 79.]& ,
       TagBoxNote->"79."],
      StyleBox["79.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[79., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{1000., 0.}, {1010., 42.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{429.0106249999999, 460.9312500000001}, {
         36.004322993146005`, 118.38926705039233`}}],
       StatusArea[#, 42.]& ,
       TagBoxNote->"42."],
      StyleBox["42.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[42., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{1010., 0.}, {1020., 30.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{459.3687500000001, 491.28937499999984`}, {
         59.09644986664496, 118.38926705039233`}}],
       StatusArea[#, 30.]& ,
       TagBoxNote->"30."],
      StyleBox["30.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[30., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{1020., 0.}, {1030., 12.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{489.72687499999984`, 521.6475}, {93.73464017689338, 
         118.38926705039233`}}],
       StatusArea[#, 12.]& ,
       TagBoxNote->"12."],
      StyleBox["12.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[12., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{1030., 0.}, {1040., 5.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{520.085, 552.0056249999998}, {107.20504751976777`, 
         118.38926705039233`}}],
       StatusArea[#, 5.]& ,
       TagBoxNote->"5."],
      StyleBox["5.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[5., {
        GrayLevel[0]}], "Tooltip"]& ]}, {}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{None, None},
  AxesOrigin->{880., 0.},
  Frame->True,
  FrameLabel->{
    FormBox["\"Steady RNA\"", TraditionalForm], 
    FormBox["\"Frequency\"", TraditionalForm]},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  PlotRange->{{All, All}, {All, All}},
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.02], 
     Scaled[0.1]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.5168769124315023`*^9, 3.5169502898125*^9, 
  3.517196243453125*^9, 3.517196514359375*^9, 3.52430462425*^9, 
  3.524398027765625*^9, 3.524902042765625*^9}]
}, Open  ]],

Cell["How many fulfil the no-drug part of the criteria :: All", "Text",
 CellChangeTimes->{{3.5167092978125*^9, 3.516709308609375*^9}, {
  3.5168769669596205`*^9, 3.516876974744066*^9}, {3.516950318375*^9, 
  3.516950320515625*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Select", "[", 
   RowBox[{"randsol", ",", 
    RowBox[{
     RowBox[{
      RowBox[{"Abs", "[", 
       RowBox[{"#", "[", 
        RowBox[{"[", "2", "]"}], "]"}], "]"}], "\[LessEqual]", 
      RowBox[{".3", "RtCheck1"}]}], "&"}]}], "]"}], "//", "Length"}]], "Input",
 CellChangeTimes->{{3.516709385875*^9, 3.51670938765625*^9}, {
   3.516710040734375*^9, 3.516710073640625*^9}, {3.51671364178125*^9, 
   3.516713643984375*^9}, {3.516714149203125*^9, 3.516714154515625*^9}, {
   3.516714190265625*^9, 3.516714209109375*^9}, {3.516876335138483*^9, 
   3.5168763826391993`*^9}, {3.5168764402084923`*^9, 
   3.5168764427106357`*^9}, {3.5168765104465103`*^9, 3.516876512173609*^9}, 
   3.516876985802699*^9}],

Cell[BoxData["1000"], "Output",
 CellChangeTimes->{3.5168764436856914`*^9, 3.5168765130936613`*^9, 
  3.516876987192778*^9, 3.5169503136875*^9, 3.524304642859375*^9, 
  3.524398035671875*^9, 3.524902048375*^9}]
}, Open  ]],

Cell["\<\
Distribution of max values of deviations from the optimal drug treatment \
kynetics:\
\>", "Text",
 CellChangeTimes->{{3.516876648613413*^9, 3.516876660539095*^9}, {
  3.516877079699069*^9, 3.5168771194983454`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"fig2", "=", 
  RowBox[{"Histogram", "[", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"Max", "[", 
       RowBox[{"Abs", "[", 
        RowBox[{"Flatten", "[", "#", "]"}], "]"}], "]"}], "&"}], "/@", 
     RowBox[{"randsol", "[", 
      RowBox[{"[", 
       RowBox[{"All", ",", "3"}], "]"}], "]"}]}], ",", 
    RowBox[{"FrameLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
      "\"\<Max. deviation from drug treatment kynetics\>\"", ",", 
       "\"\<Frequency\>\""}], "}"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.516876664801338*^9, 3.5168766905718126`*^9}, {
  3.5168771846420717`*^9, 3.516877192743535*^9}, {3.516877234937948*^9, 
  3.5168772377591095`*^9}, {3.5168772739081774`*^9, 3.516877304671937*^9}, {
  3.517196281765625*^9, 3.51719631271875*^9}, {3.517196519*^9, 
  3.517196519953125*^9}, {3.524387316109375*^9, 3.52438731890625*^9}}],

Cell[BoxData[
 GraphicsBox[
  {RGBColor[0.798413061722744, 0.824719615472648, 0.968322270542458], 
   EdgeForm[Opacity[0.7]], {}, 
   {RGBColor[0.798413061722744, 0.824719615472648, 0.968322270542458], 
    EdgeForm[Opacity[0.7]], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{0., 0.}, {0.5, 16.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{64.71312499999999, 80.56180147058822}, {
         106.00883374028471`, 118.38926705039233`}}],
       StatusArea[#, 16.]& ,
       TagBoxNote->"16."],
      StyleBox["16.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[16., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{0.5, 0.}, {1., 17.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{78.99930147058822, 94.84797794117645}, {
         105.332712908403, 118.38926705039233`}}],
       StatusArea[#, 17.]& ,
       TagBoxNote->"17."],
      StyleBox["17.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[17., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{1., 0.}, {1.5, 34.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{93.28547794117645, 109.13415441176468`}, {
         93.83865876641366, 118.38926705039233`}}],
       StatusArea[#, 34.]& ,
       TagBoxNote->"34."],
      StyleBox["34.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[34., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{1.5, 0.}, {2., 12.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{107.57165441176468`, 123.42033088235291`}, {
         108.71331706781163`, 118.38926705039233`}}],
       StatusArea[#, 12.]& ,
       TagBoxNote->"12."],
      StyleBox["12.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[12., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{2., 0.}, {2.5, 7.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{121.85783088235291`, 137.70650735294114`}, {
         112.09392122722025`, 118.38926705039233`}}],
       StatusArea[#, 7.]& ,
       TagBoxNote->"7."],
      StyleBox["7.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[7., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{2.5, 0.}, {3., 13.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{136.14400735294114`, 151.99268382352938`}, {
         108.0371962359299, 118.38926705039233`}}],
       StatusArea[#, 13.]& ,
       TagBoxNote->"13."],
      StyleBox["13.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[13., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{3., 0.}, {3.5, 10.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{150.43018382352938`, 166.2788602941176}, {
         110.06555873157507`, 118.38926705039233`}}],
       StatusArea[#, 10.]& ,
       TagBoxNote->"10."],
      StyleBox["10.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[10., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{3.5, 0.}, {4., 8.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{164.7163602941176, 180.56503676470584`}, {
         111.41780039533853`, 118.38926705039233`}}],
       StatusArea[#, 8.]& ,
       TagBoxNote->"8."],
      StyleBox["8.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[8., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{4., 0.}, {4.5, 16.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{179.00253676470584`, 194.85121323529407`}, {
         106.00883374028471`, 118.38926705039233`}}],
       StatusArea[#, 16.]& ,
       TagBoxNote->"16."],
      StyleBox["16.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[16., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{4.5, 0.}, {5., 12.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{193.28871323529407`, 209.1373897058823}, {
         108.71331706781163`, 118.38926705039233`}}],
       StatusArea[#, 12.]& ,
       TagBoxNote->"12."],
      StyleBox["12.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[12., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{5., 0.}, {5.5, 24.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{207.5748897058823, 223.42356617647053`}, {
         100.59986708523091`, 118.38926705039233`}}],
       StatusArea[#, 24.]& ,
       TagBoxNote->"24."],
      StyleBox["24.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[24., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{5.5, 0.}, {6., 22.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{221.86106617647053`, 237.70974264705876`}, {
         101.95210874899436`, 118.38926705039233`}}],
       StatusArea[#, 22.]& ,
       TagBoxNote->"22."],
      StyleBox["22.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[22., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{6., 0.}, {6.5, 32.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{236.14724264705876`, 251.995919117647}, {
         95.1909004301771, 118.38926705039233`}}],
       StatusArea[#, 32.]& ,
       TagBoxNote->"32."],
      StyleBox["32.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[32., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{6.5, 0.}, {7., 42.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{250.433419117647, 266.2820955882352}, {
         88.42969211135984, 118.38926705039233`}}],
       StatusArea[#, 42.]& ,
       TagBoxNote->"42."],
      StyleBox["42.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[42., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{7., 0.}, {7.5, 79.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{264.7195955882352, 280.5682720588235}, {
         63.41322133173598, 118.38926705039233`}}],
       StatusArea[#, 79.]& ,
       TagBoxNote->"79."],
      StyleBox["79.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[79., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{7.5, 0.}, {8., 237.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{279.0057720588235, 
         294.8544485294117}, {-43.413870105576734`, 118.38926705039233`}}],
       StatusArea[#, 237.]& ,
       TagBoxNote->"237."],
      StyleBox["237.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[237., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{8., 0.}, {8.5, 407.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{293.2919485294117, 
         309.1406249999999}, {-158.35441152547014`, 118.38926705039233`}}],
       StatusArea[#, 407.]& ,
       TagBoxNote->"407."],
      StyleBox["407.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[407., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{8.5, 0.}, {9., 5.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{307.5781249999999, 323.42680147058815`}, {
         113.4461628909837, 118.38926705039233`}}],
       StatusArea[#, 5.]& ,
       TagBoxNote->"5."],
      StyleBox["5.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[5., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{9., 0.}, {9.5, 2.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{321.86430147058815`, 337.7129779411764}, {
         115.47452538662888`, 118.38926705039233`}}],
       StatusArea[#, 2.]& ,
       TagBoxNote->"2."],
      StyleBox["2.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[2., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{9.5, 0.}, {10., 1.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{336.1504779411764, 351.9991544117646}, {
         116.1506462185106, 118.38926705039233`}}],
       StatusArea[#, 1.]& ,
       TagBoxNote->"1."],
      StyleBox["1.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[1., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{10., 0.}, {10.5, 1.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{350.4366544117646, 366.2853308823528}, {
         116.1506462185106, 118.38926705039233`}}],
       StatusArea[#, 1.]& ,
       TagBoxNote->"1."],
      StyleBox["1.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[1., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{10.5, 0.}, {11., 1.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{364.7228308823528, 380.57150735294107`}, {
         116.1506462185106, 118.38926705039233`}}],
       StatusArea[#, 1.]& ,
       TagBoxNote->"1."],
      StyleBox["1.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[1., {
        GrayLevel[0]}], "Tooltip"]& ], {}, 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{11.5, 0.}, {12., 1.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{393.2951838235293, 409.14386029411753`}, {
         116.1506462185106, 118.38926705039233`}}],
       StatusArea[#, 1.]& ,
       TagBoxNote->"1."],
      StyleBox["1.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[1., {
        GrayLevel[0]}], "Tooltip"]& ], {}, {}, {}, {}, {}, {}, {}, {}, {}, 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{16.5, 0.}, {17., 1.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{536.1569485294116, 552.0056249999998}, {
         116.1506462185106, 118.38926705039233`}}],
       StatusArea[#, 1.]& ,
       TagBoxNote->"1."],
      StyleBox["1.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[1., {
        GrayLevel[0]}], "Tooltip"]& ]}, {}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{None, None},
  AxesOrigin->{0., 0},
  Frame->True,
  FrameLabel->{
    FormBox[
    "\"Max. deviation from drug treatment kynetics\"", TraditionalForm], 
    FormBox["\"Frequency\"", TraditionalForm]},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  PlotRange->{{All, All}, {All, All}},
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.02], 
     Scaled[0.1]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.5168769124315023`*^9, 3.5168773176386786`*^9, 
  3.5169503310625*^9, 3.51719632371875*^9, 3.5171965241875*^9, 
  3.524304678046875*^9, 3.524383670421875*^9, 3.52439805471875*^9, 
  3.52490205575*^9}]
}, Open  ]],

Cell["How many fulfil the drug part of the criteria:", "Text",
 CellChangeTimes->{{3.5167092978125*^9, 3.516709308609375*^9}, {
  3.5168769669596205`*^9, 3.516876974744066*^9}, {3.516877170613269*^9, 
  3.51687717097429*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Select", "[", 
   RowBox[{"randsol", ",", 
    RowBox[{
     RowBox[{
      RowBox[{"Max", "[", 
       RowBox[{"Abs", "[", 
        RowBox[{"Flatten", "[", 
         RowBox[{"#", "[", 
          RowBox[{"[", "3", "]"}], "]"}], "]"}], "]"}], "]"}], "\[LessEqual]",
       ".27"}], "&"}]}], "]"}], "//", "Length"}]], "Input",
 CellChangeTimes->{{3.516709385875*^9, 3.51670938765625*^9}, {
   3.516710040734375*^9, 3.516710073640625*^9}, {3.51671364178125*^9, 
   3.516713643984375*^9}, {3.516714149203125*^9, 3.516714154515625*^9}, {
   3.516714190265625*^9, 3.516714209109375*^9}, {3.516876335138483*^9, 
   3.5168763826391993`*^9}, {3.5168764402084923`*^9, 
   3.5168764427106357`*^9}, {3.5168765104465103`*^9, 3.516876512173609*^9}, 
   3.516876985802699*^9, {3.5168771576775293`*^9, 3.516877163702874*^9}}],

Cell[BoxData["2"], "Output",
 CellChangeTimes->{3.5168764436856914`*^9, 3.5168765130936613`*^9, 
  3.516876987192778*^9, 3.5168771648989425`*^9, 3.516950346890625*^9, 
  3.52438368075*^9, 3.524398065828125*^9, 3.5249020675625*^9}]
}, Open  ]],

Cell["Check only two drug concentrations:", "Text",
 CellChangeTimes->{{3.524387161203125*^9, 3.5243871774375*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"fig2", "=", 
  RowBox[{"Histogram", "[", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"Max", "[", 
       RowBox[{"Abs", "[", 
        RowBox[{"Flatten", "[", "#", "]"}], "]"}], "]"}], "&"}], "/@", 
     RowBox[{"randsol", "[", 
      RowBox[{"[", 
       RowBox[{"All", ",", "3", ",", 
        RowBox[{"{", 
         RowBox[{"1", ",", "2"}], "}"}]}], "]"}], "]"}]}], ",", 
    RowBox[{"FrameLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
      "\"\<Max. deviation from drug treatment kynetics\>\"", ",", 
       "\"\<Frequency\>\""}], "}"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.516876664801338*^9, 3.5168766905718126`*^9}, {
  3.5168771846420717`*^9, 3.516877192743535*^9}, {3.516877234937948*^9, 
  3.5168772377591095`*^9}, {3.5168772739081774`*^9, 3.516877304671937*^9}, {
  3.517196281765625*^9, 3.51719631271875*^9}, {3.517196519*^9, 
  3.517196519953125*^9}, {3.52438721190625*^9, 3.524387219765625*^9}, {
  3.524387324984375*^9, 3.5243873279375*^9}}],

Cell[BoxData[
 GraphicsBox[
  {RGBColor[0.798413061722744, 0.824719615472648, 0.968322270542458], 
   EdgeForm[Opacity[0.7]], {}, 
   {RGBColor[0.798413061722744, 0.824719615472648, 0.968322270542458], 
    EdgeForm[Opacity[0.7]], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{0., 0.}, {0.5, 16.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{64.71312499999998, 80.99471590909089}, {
         108.34333957984549`, 118.38926705039232`}}],
       StatusArea[#, 16.]& ,
       TagBoxNote->"16."],
      StyleBox["16.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[16., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{0.5, 0.}, {1., 19.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{79.43221590909089, 95.7138068181818}, {
         106.75269692911797`, 118.38926705039232`}}],
       StatusArea[#, 19.]& ,
       TagBoxNote->"19."],
      StyleBox["19.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[19., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{1., 0.}, {1.5, 32.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{94.1513068181818, 110.43289772727269`}, {
         99.85991210929868, 118.38926705039232`}}],
       StatusArea[#, 32.]& ,
       TagBoxNote->"32."],
      StyleBox["32.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[32., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{1.5, 0.}, {2., 14.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{108.87039772727269`, 125.1519886363636}, {
         109.40376801366385`, 118.38926705039232`}}],
       StatusArea[#, 14.]& ,
       TagBoxNote->"14."],
      StyleBox["14.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[14., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{2., 0.}, {2.5, 10.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{123.5894886363636, 139.8710795454545}, {
         111.52462488130055`, 118.38926705039232`}}],
       StatusArea[#, 10.]& ,
       TagBoxNote->"10."],
      StyleBox["10.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[10., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{2.5, 0.}, {3., 10.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{138.3085795454545, 154.59017045454542`}, {
         111.52462488130055`, 118.38926705039232`}}],
       StatusArea[#, 10.]& ,
       TagBoxNote->"10."],
      StyleBox["10.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[10., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{3., 0.}, {3.5, 12.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{153.02767045454542`, 169.30926136363632`}, {
         110.46419644748221`, 118.38926705039232`}}],
       StatusArea[#, 12.]& ,
       TagBoxNote->"12."],
      StyleBox["12.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[12., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{3.5, 0.}, {4., 10.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{167.74676136363632`, 184.02835227272723`}, {
         111.52462488130055`, 118.38926705039232`}}],
       StatusArea[#, 10.]& ,
       TagBoxNote->"10."],
      StyleBox["10.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[10., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{4., 0.}, {4.5, 13.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{182.46585227272723`, 198.74744318181814`}, {
         109.93398223057302`, 118.38926705039232`}}],
       StatusArea[#, 13.]& ,
       TagBoxNote->"13."],
      StyleBox["13.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[13., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{4.5, 0.}, {5., 16.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{197.18494318181814`, 213.46653409090902`}, {
         108.34333957984549`, 118.38926705039232`}}],
       StatusArea[#, 16.]& ,
       TagBoxNote->"16."],
      StyleBox["16.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[16., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{5., 0.}, {5.5, 26.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{211.90403409090902`, 228.18562499999996`}, {
         103.04119741075374`, 118.38926705039232`}}],
       StatusArea[#, 26.]& ,
       TagBoxNote->"26."],
      StyleBox["26.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[26., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{5.5, 0.}, {6., 26.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{226.62312499999996`, 242.90471590909084`}, {
         103.04119741075374`, 118.38926705039232`}}],
       StatusArea[#, 26.]& ,
       TagBoxNote->"26."],
      StyleBox["26.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[26., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{6., 0.}, {6.5, 42.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{241.34221590909084`, 257.6238068181818}, {
         94.55776994020691, 118.38926705039232`}}],
       StatusArea[#, 42.]& ,
       TagBoxNote->"42."],
      StyleBox["42.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[42., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{6.5, 0.}, {7., 50.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{256.0613068181818, 272.34289772727266`}, {
         90.3160562049335, 118.38926705039232`}}],
       StatusArea[#, 50.]& ,
       TagBoxNote->"50."],
      StyleBox["50.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[50., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{7., 0.}, {7.5, 130.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{270.78039772727266`, 287.0619886363636}, {
         47.89891885219939, 118.38926705039232`}}],
       StatusArea[#, 130.]& ,
       TagBoxNote->"130."],
      StyleBox["130.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[130., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{7.5, 0.}, {8., 519.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{285.4994886363636, 
         301.7810795454545}, {-158.3544115254702, 118.38926705039233`}}],
       StatusArea[#, 519.]& ,
       TagBoxNote->"519."],
      StyleBox["519.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[519., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{8., 0.}, {8.5, 49.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{300.2185795454545, 316.5001704545454}, {
         90.84627042184268, 118.38926705039232`}}],
       StatusArea[#, 49.]& ,
       TagBoxNote->"49."],
      StyleBox["49.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[49., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{8.5, 0.}, {9., 2.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{314.9376704545454, 331.2192613636363}, {
         115.76633861657396`, 118.38926705039232`}}],
       StatusArea[#, 2.]& ,
       TagBoxNote->"2."],
      StyleBox["2.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[2., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{9., 0.}, {9.5, 2.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{329.6567613636363, 345.93835227272723`}, {
         115.76633861657396`, 118.38926705039232`}}],
       StatusArea[#, 2.]& ,
       TagBoxNote->"2."],
      StyleBox["2.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[2., {
        GrayLevel[0]}], "Tooltip"]& ], 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{9.5, 0.}, {10., 1.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{344.37585227272723`, 360.6574431818181}, {
         116.29655283348315`, 118.38926705039232`}}],
       StatusArea[#, 1.]& ,
       TagBoxNote->"1."],
      StyleBox["1.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[1., {
        GrayLevel[0]}], 
      "Tooltip"]& ], {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, 
    TagBox[
     TooltipBox[
      TagBox[
       DynamicBox[{
         FEPrivate`If[
          CurrentValue["MouseOver"], 
          EdgeForm[{
            GrayLevel[0.5], 
            AbsoluteThickness[1.5], 
            Opacity[0.66]}], {}, {}], 
         RectangleBox[{16., 0.}, {16.5, 1.}, "RoundingRadius" -> 0]},
        ImageSizeCache->{{535.724034090909, 552.0056249999999}, {
         116.29655283348315`, 118.38926705039232`}}],
       StatusArea[#, 1.]& ,
       TagBoxNote->"1."],
      StyleBox["1.`", {
        GrayLevel[0]}, StripOnInput -> False]],
     Annotation[#, 
      Style[1., {
        GrayLevel[0]}], "Tooltip"]& ]}, {}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{None, None},
  AxesOrigin->{0., 0},
  Frame->True,
  FrameLabel->{
    FormBox[
    "\"Max. deviation from drug treatment kynetics\"", TraditionalForm], 
    FormBox["\"Frequency\"", TraditionalForm]},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  PlotRange->{{All, All}, {All, All}},
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.02], 
     Scaled[0.1]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.5243873361875*^9, 3.524398070390625*^9, 
  3.524902072234375*^9}]
}, Open  ]],

Cell["\<\
The no-drug steady RNA vs. the maximal deviation for the drug treatment:\
\>", "Text",
 CellChangeTimes->{{3.516877617764845*^9, 3.5168776453494225`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"fig3", "=", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"ListPlot", "[", 
      RowBox[{
       RowBox[{"MapThread", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"{", 
           RowBox[{"#1", ",", "#2"}], "}"}], "&"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"RtCheck1", "-", 
            RowBox[{"randsol", "[", 
             RowBox[{"[", 
              RowBox[{"All", ",", "2"}], "]"}], "]"}]}], ",", 
           RowBox[{
            RowBox[{
             RowBox[{"Max", "[", 
              RowBox[{"Abs", "[", 
               RowBox[{"Flatten", "[", "#", "]"}], "]"}], "]"}], "&"}], "/@", 
            RowBox[{"randsol", "[", 
             RowBox[{"[", 
              RowBox[{"All", ",", "3"}], "]"}], "]"}]}]}], "}"}]}], "]"}], 
       ",", 
       RowBox[{"FrameLabel", "\[Rule]", 
        RowBox[{"{", 
         RowBox[{"\"\<Steady RNA\>\"", ",", "\"\<Max deviation\>\""}], 
         "}"}]}]}], "]"}], ",", "\[IndentingNewLine]", 
     RowBox[{"Graphics", "[", 
      RowBox[{"{", 
       RowBox[{"Red", ",", 
        RowBox[{"Line", "[", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"800", ",", ".27"}], "}"}], ",", 
           RowBox[{"{", 
            RowBox[{"1200", ",", ".27"}], "}"}]}], "}"}], "]"}]}], "}"}], 
      "]"}]}], "}"}], "//", "Show"}]}]], "Input",
 CellChangeTimes->{{3.516877652340822*^9, 3.516877772450692*^9}, {
  3.516877806358632*^9, 3.5168778637719154`*^9}, {3.516950363875*^9, 
  3.516950368125*^9}, {3.51719653290625*^9, 3.517196533921875*^9}}],

Cell[BoxData[
 GraphicsBox[{{{}, 
    {Hue[0.67, 0.6, 0.6], PointBox[CompressedData["
1:eJwVmnk4VW0UxYWQKWPGe88xVEhSoUzdTaaQKFSSSpJERVEiRGW651xERZEi
FJGpQUREElJJkqEoJEmSJPTt76+e+9yc4d17rfVbD0ruhzfv4+bi4hri4eL6
/9/L4UvqrjVzQDt3bZudbCtr3rp+JX+ahuURD71aAxiQXzFd4PiJBqeyaOPM
r5JQ565CDS6m4VerjsXtYRlw27POpqGBAgn+q++3ppAQ9172pJUpDRe2mNWP
3lCGymvWBwo6adDfI/B33hMCilpO+IQlUbB/1qEsejsB5Yz3MbrlNKhEpMF1
SznwN1TbrHCbgrVWj5cYiP1jyYZsaFo0TMOWj2fXcHErwO4+gQLCmwYFXu79
a+8SMH5wbCKIiwN8ynJPE9eTEOGh96tEnwMLDE5FDfwmoH2yJum4Qjxkffo2
3HdwhNXRHK5atIWCT3pyaafdSLgv9Cr4I77PTNDWU2sd8OfHDCI9Myi499VD
at12Rdh8ruBkVzQFeYHxCqNLpCDPeyHkeFCQUmZofcJcATQjKY/P5Rz4Vftl
yEWKAFr4pnKhJg2woYiSAxLin0tsdvnKgf4vU9yrfjEgsMXl6gEHGn7bg0OU
nhxoNP+4dsmMhr6ulIttYUvhlr3Yd/7bNLTNHdvqfHGCpdTBbHt1jIIAozsL
RAwIsO4Tvx2eSoPR2UVHNpmQ4GP3VuKzEw13g4iADnz/sS6pK9G2NAySqtIJ
c4rgwjYVWJdNQd+LveuX8ZMwed31qX8dDVEiVWMVxYsg58KsgGQpBRk/JkpZ
wj2s/nPXlY6+oUD4mMiLwEX8oF26vDjDhQbNE1x3ogVJ+PSztvJAEQfuBS5I
Fr6iCFcorqhMfQrEttWoTnqQIGgUfC52nIIg36c/r71kglR+g9Y8moJd5YUt
W9/JgGfhi01+/8/nX9mA6i5F2G7k0A38NHRc5qxWPEhC5o97E5I4n4un9707
V6YALU2zxZ7hNFjs0AkdcSRheqRP9sg8GqyWaaz0O0+C6dLk5Bg3DvyZVGp0
DiTAqvKOZP03CtYNFAp6v51k6W2npql/FDQpuhrXmP9jWboHhXot5QCcexqT
uBTnrZfCc7yPDVI21WmJR0iw5q+dSxWjYce8+F+ru5lwQZXHkruOAq+ejYqK
aiQ8q/FjWLfg/S+mDueofGbxDGnvjfvAhrsxKXoMJxJuT58+rd1JwaSLbvXc
6Djr0OOH+/rnx4Pwwb7K2g1CwLXwh7PJZRrU6lxfjZ0lIe3iC3fxQQp0Xx7O
nLIQg1Wlaq/24747CWzo5+tkQsXBFdPL1tLQvLsmQOcsE4KcmBL7xChY+bTH
/PktJqwaMnUKyaOh18R/7WXZRTB0J+N4uhcNQWmMPXBRDmx3RZULy1Ew/+oz
96AVqLcEJ4vHFjRovP96zu00CUf8/4zcaKKgzAYcubQJaOrTWp1/kQLBAJnw
YBESbr2SC+Vm0fBXbZvh3EoGaArvmjOppyErov7Gtp0kxDY3qGZbcuCTz5Xu
hRtJOBRT2Bi0n4IuWa7CDNzP/H/vrz+6woEgYXJavIYA+5MK39joDx9fOP26
e0wKjPyfec/NceChaY/xZLEiCCYOD3/4QMHQjfqy1CIShIw/xP02okFdQKZh
yRFFeOIXklM6SwH3nctCCqjfYxN3hjZmsoH84lpZivuzXujAC4UEDtTZTZdX
vmBCdVNe1ewCDlQwmB3NW36x+hI7u5bdo2Eg9OhgpzkJ/rszpR69Z8OO89mv
BA0J2Cab8SIB933Xis27U44QELLvAtd+Ew6wzXMbVVHPIsYG+t6NeH6X/l7K
ElGCrrZl4WGxHPB9U1/k4s2EYTEFvRYmBxYTpQvXpBIw8a0mRvMdDSYb9l5V
qmTCnZZwh7vdNJQ179gUkkTA5+Ka04e8KVj8W3TvL2NpiPjowjuzj4Kry37v
FOGThj7g3M8/TsPq0nklYqvk4LKS68oqPL9rDKlnoSFMKB/a9HK1PwUnd1X7
bBFYCI+C5k7urqSBdTk8Rr+WAW7EMh4X1M+RCbWPO4JqWR6SiSVFgTS0L1vM
lbyIhKHtA2TGJBvcajj0ZisSWqTk/+TksIHvnFXSNwMSws72/XU0ZEOQq5ey
bcl8mAmpv2HRTsF9j5BNS8NJgGi/gVqg4ao8Ta3cw4CQ+uupPPi+D38p2syl
ETDs9sonXpMC8zWWE01vCWC7jVdlOrChW9c0LugrA5p2J+rqbeZAnt/HaOKv
DITyjAiUnaUgfiS5TD+vm+VYw+hz1kO9Pnk2eOZqD+vW5/4ocW0OFAR2yivL
CcMg9+D+oEWo94sVGqp2JFxNow/0KNKw9cEXedUzilAg93b9BksK3MQquWeX
kKCy0Jhr+S08P96IFCMTAq7uYOT6a2BeRV3f8x3ne9vXN9coiwNix7sPnkI9
c4mQym/SKSgvGiy08Sfh8MarjtFtNJyeWP99pEQGHi3QVRBspWBf5Jtze5sU
QbRDymzuPO7n7ZVqA2YkGNouu9D2F99/2P8S3UVAhFiMrFYMziPmp7/VJxKY
/CJ+y77TML2ccVEN/UhqZbN6F/r/KWcbdytvAvY+FlioeBzz5+yVGAH087gD
zpZXP7OhS+fvm7d4/qFPgp38drIBrr4rPYLPx+znts7w4oCawFD+fNcpVpzj
uaOiJTRsPDZ2QFqLAOEv6itXhtFwftlU8eMLiyB0s33usnQadm4T4G9XVYSK
JVlJk01sqLikP6uN1/u0rWCly0UaYo1Lvw9uwrxJtA1fxkfD87PeN6JmcJ5y
eg1dJRSIijSdaj5Gwv4XptneZngeA+8UojQIKNVjsTUxT3/F/vLoipKA5e9F
jVLwvBaoq9XXNhCwtsiG504/BVorbWePKZGQnuDCl66F/ugRvWSW4gOJLclv
uzHvFIau7W1JIMD92B0vcysKbj/lHAhZJgn5GVnD68Q48LizwuCUJhOexDdN
Xexnw6+UCW2ePwT4r4j1iiyngOdpY+Z+GQLEVGpjU11pmDpntP6VNgnXD7w5
MH+Whqbs75Vnz31hXW9oeRKRSUFkgvuGW2wCUus2vvRQpCDhsc/skauKoElt
Hp/dTUNab+K/L3lMMPSrN6wx5oBVlPOV+cok6BVe6Srz5MCNiqW1pz5KA6kd
ae7+i4Ie3S8au8VJ6IuveFONeS8sJVN4JoMEAaHSm1QXDVVewd/5ZplAqxz/
fusjBXdiJlSPaJEQ5F4gNv89Ddo2DzJHF5Ow9uES3SDcx+AO36zL4/Jw6Gns
pvXoB9aWS8vGFEio+PajcnqAAslMrchDzgRMz765cMAR/XYs8eN8Llmofnyo
90MGDVfq91NxW0ioHz9U/wzzaLmIcif9iYCZf3cMLZwp0Gm513RmD+bjniiF
2RAaXt1vvlGEz7M+/rk0z0saMm5Ntq+SI8FMJkN0vyMFe9V+tJbVElBMeW86
qEfBzgUP+n20FCBDObzMYpyGn4LrT3x7KgE3N44EGIvQMGGnVJCHedHA0xg9
Is+BtQdtnpxCP1q363SdiSE+z+eHNdVH5KD4+omLCY8p8Bddc8PeiATW4A1Z
2514fYsgV5UsRYgi2kf0UE+/9lJbby4Sg9E/KdTOvRwourHH1r4J/dt8yI0V
TINuHfNh1kVFqF4AU1uRTxxlf2X16pDQqvnWtWoNB7RcHmTvO0RCXt6Fdz01
NHzLv/lo30MCvC6rL5xXyAaj5Oq8vsMk8HhGDQbh/Fpd6gxE8bzX6fbWl6+k
oF/rZkRHIANODC1sf4h8Z2fmYLwwhITn8nlb/jRTQDpcPzaoQULpOd0wvs8U
tOxpZJ9cjZ/pLEeFVTSc/dpUv66AgKhdtU3GKyi41DbbK7oW9eD7fJjCPF6b
993x7QrM2630PBZJgZzNvQMfKpiw8kCa8sbtFLCnrf5s207CwtOLnfwTaTjw
6M7ZCJoB3iLhwHxCgd225aVLSRKuMPrWNAWjPsyi6VEVWTjRJu9TFYu8cW79
VPUWJnQkO3+29KOhyHi5efwvWbDV3vFSQZIC5637FpVyk+A1eIu0aaeBT19/
42fMi1uLE4Zkr1Jw4p2g4IfzI6x9Q3VWwV/i4J73UZVVC2RgufaFIDv0f68X
hYHPri0A+cxwax/kbWbkle8uqL9DAefP1mDee5T2igmtI2FP4O8PzzwpEGkZ
bjNYiPlZ4HDdrQ1529dBvtWRCRtE84L1NtDQcE72LO1OwpTRzc4N9yngX75w
xbOlCnDI45V2yVfUl+cmXn2CBKXnwlJjkeh/t9g3f24igJfnX/MY8nJF8fTF
/q0kdFn8wybCAZU/R4ZcPkkCr1GXBb8yDb6S+yySxggwihnWl7xEQX7tFsda
5OfdqSoTn2RpyLT67s7LYUIqeSaMQH/61+JNRvIrwieNwh+ROM+47waLOi/J
wrcdNk9OXsZ9cz45NrMD+fX9il1LkPfMhj2PCXmT4Ba6uET3ARsKjAPjPquQ
sGSE1SSaSUOem0WSXjYTWsMfUzd6KPiRkG2aly0PJ/uCXPOzaBCZevzt2G55
kMyb8nq8n4bPLxOk0/oJGGO4h5pZ01BXnXc3Y58s5LvKRKzi40BN48d7sahf
x7tRmklrKFj+rjF3UocJm32fdv38y4H8HUopHriP5JfCWv1UDqifdhodx341
2H/GZ2E3BQy7uGxmOgHZb9czugooMGCmCx/CfVJ86HBuyIEDO1Sf23Cy5KFa
f87zsibmZ6v5xm/VAvB7Ja/b6jQOUG0JZTy5ipD9c2K33HIa5soVjksvJ0Fz
p1TgsQU0vODatPHuJAHHB08O9HBzQOZMI2/NSgkQ/JbX7Yr8nfJ7Q5v0azEQ
bqdc7dvZYJi+/LmdOwF8Nu9nWt+w4X51T6nXCRJOenIPJKFfWLfoT+tak9De
UHC9HPlvQmpdTosNCWJMiz9f5NlQWhLhGoZ9L7Xt+aL+Fgq2dqaWiyIvSknd
FVu7lAZS4KbtU+Sf5CM3y/aexzwWiudqvioPx+YZZUZP0XDoyKaGk+pjLMH9
H59rSXHA1LVQ5XXPfIiytpUzwXyNlza6EBAzwspvNovZdogNcV49u4oCmKBd
IXqRGqFAuvRopt0UAScv5sSdmUAedOxylz1JwK1Yh68rHtHAq5Rbkr4K/Tb1
4dcTEjR4/plu00YeEc4+XHD9N/I68ZsaucEAwaONam9zadjDnogRLCLg0Ido
F11/3Id3Y/wZywlIawgRSEKebFlu23jSRBgK1A+HHduHfqPbmGSZSCC3bpGy
FkR/PPHl4jL0/wq+jq68v2yoFM091ZBBQPydxSTPRzZYbLasP4Ln4Tl5xqXe
GP3k2hv7D8jPfB9rZ01xnlG5lUIDWQw4s1RrcUgsBXtSPhbxxDOhxL6rVfQj
DZPndqZq65PwsU1sgAfPe/kihwduWwgwyP1pdnYqDshkboU1sUJgx2yNYk1R
0NkYwJ9ujPt6qUZiO/rnit5THQIFDPBnDqamIB9ZLPELcsF5zxQflsoOxv2Q
9s+xwPwvfbj0dST68clHCSu5MS9dyi7uNET/YLKy+j3Rb94y7W0cOjjwetDy
GleWGBz0SeskDCjQ19FNmKphQs3T6EO68TQscDncfWGOAczbd798W47+emR/
sbCGEsBUkE0M7peKo32/T6kidDTdc37mgH135NJIjxzyxsSm65eQ96r2KD19
ZioOXNGP32vSNEg+DPRyP8OEvUn2ri0qHJh/TTX9nLA80Ke2HlFEntixr9Hn
vi/u897u9s18FIgPDp0xQL2NZ3Vf7drIhpnaJYlTC0kwyNo5+5GbhgQ1I7+z
2B+cpLqt7JNo8Hokc9SufTFIdNavd/pAQ99tvZzcOFX4NCAnNFBJgevQo8M6
mJ+d654vb1qK/a7x9RMVQ+R965Cc+FEaIqJfmmQif8ZV+RF/kFeiOB6p5wYZ
kL1MfYf/CA31BaJ30vD+DceGqxecoMGPI/syV0gOro/+kfDdS8H3z331p4EJ
EWtT9tzB/PdLiu6uPKwAgRqbKK1s3G/tXb6tSUy4/v3uU4lezAOuS4LeVoKg
0Vf3668hB454Z7p0rCFBq65z11peGu5k5jx/E9TDEo1qTriJvJMisHBbbDED
rt+OLZVzZUN5MqXFwX6pKNIcuhL9LWSm6+Sb/SSMD6R570O/MRtRv1GxdohV
kHhiw4pv2E92E661fXIQH/hZQAf1LV7SwUVRsiCqr/0zAffRk7g1HJmK/V7P
L9dXHftjdeGahzoy4HbWh/3gPg0X3RKT2nzkgd9Hu+uiOwVN/4jQewKfWPVe
BV/ebKZA81KinZz2JEt5TWEqGY36NYpRuo99wilQprAP+3PZcr9ayS8SIFQi
ZWelTUFpiHXZIlUSkoq38+09RoPqusQt4cg/A2O9G1agXlL81gXs7JGDtofC
6v2qmEfj7aM/jsyDASepTa6Yp9mbmVCFfeb+g+8ZsrcpmNaSOfMU+/AOgXBL
QexHjK3hVecFCVhEaJncw7xV/sm1vwr1V36ze/E08t66+/w2KhYkVP/r3Z0u
RYNW2HYdf8wjj4fcmunYfxWjnUVcMI8PxApGfEqjwYznxDtv5DW1++0O+xNo
oBM3xu+TIUH69nnJJXE0jBfXMM8jn9j5F6o3FVIQRzxSMXNjgPS4R9G/n7hP
H35keaPfXfn257HvfexvLrXzcs+SYLR1cV6VCPbrhxPlUshHF1Rep54Y5cDZ
BWPCDUwm/LnxWzIW80Cz2mdfEvLmsUMbDcwmKchS+9PLEkN9X2uiW5DHakfX
1l9HfhJ0SeJvKqIgObBLh/OeANrFPYYL581tu/idH4eEpVs1HmgZsMG+cS56
myIJonsPBr+w54Br5pvYCUfMYx85q0rct0ema2985UV9pA5WXsY+2XDrX0fX
PuTlofR7R1+wgTSf9/CcKwnfHP+pzzvGBsJYJ8Dugjh45F+nLNG/ZBZ9yzwd
JA5GZLF5/DANBvsHr53F8/eYF+hajP148+NyoQt4Xkvrn/SvkabBXfzsDX1T
fD82PTQQSIHvn7ivdXvnQ5ecUPL0QwqijZY5GduSoDNcw1WCPDaZfPvPR8kf
LEbsZGMB8tLNqHDutELkJfpJbCX6Uc7vkcLh64oQ0RJjNeFMw/0rweXRXMjz
xw9vX4/7m5GzbLoM/fGXipBTPeZH9fFixTYtZXgWb2G14Q8Htswue24rxoTY
4ky3bOSP9DSlr7StAAhPn+b0MSh4U9Cw7tQqbviqUiP5CflrW39Oc3JnD0s3
stftzTjyZb6d90vkN1PRUIO60zT0BumFqpTJQF/iytfnDtIwyrHT+aOJ885Q
zR6lsW9YHdUXw/d34LpjlYp+9tLn5tCJbUyoPvxm65WjqM80z2vb0R80PnPv
uCJDw94Sh7dWIthn14a93G1Ag+tX9V28uE+vpc8tTO2mQWNBn50J8nCI22Ay
Vx8NhZNTE9cZ8nB39j73FpznmfwSdybqxch525loPK8Tl7RuXHSaZLn/dD5a
00VB6gAz9oMg8vvppwqjyDNvnCwYVJE03Dl946HSVhp2h18Pvo76GNPYvL0r
hwPZJo783chXrEy/DTtQz6VLOa8bhcdZHpai+ld+01BebF6m9/of6+69D8Ny
nRS8dope/NlCEqI3nL+bgvyjYv9e4+1GBdCzD1je9owNqhp3H+exSbi7+hpj
APuw5O6dDn69TMirvc0vZEGDnECC32pzOQiua2/o3Yj9JzotUicG+2r8zpKf
FRQMuzeZZpzB87owSXeMYZ9pWJf77wsTVqw4RARhP7//VljzO/qn66exSn3k
KesrUq4pH6dYggd031V/oeCGSkSN/FJp4Gcb3tiP/rZ2/q9gJ/SLxNAbC8p/
xEEjJ/MsrUdCkYOYprkRDUcOanzR8kQ/eygX1orzVeD4bVyBfPGS672g4AwF
oYE+xfYHSCjeKk4+OkxDz6tn13VzMW/CHY6xCBq+8Gra2rfgPLUPJCff4sB6
XmHyQgoD4gOC3n+5jv019pt52E9FGOs9VJC4iYKzJ2IWnHjMhPFz8Xue3KXA
/c5vez18vp3Kk+/DkH8U3rxb07iOAb0lY0RAPxu2LdCfDcV9WLZhhb2pWjyQ
3n2H0jdh3gh43wpK5wCZstpjenCcZd15yng9zuPHhDmfZ4EMeLzcmuOPfdBY
8uz5XuxTS45V+1Q9o+Bi8j4R6ZuY19IN+TcD2KCbp33apYeAHx9ntficKHhc
bhywELjBTkdrrsKMAx1O1xOVhWdYlX6qny9Mo97MeJM6OiZYk9/7n7tr0sC/
Y+5XrfYiGBUaT+D6gnkTHKGSPEOARcQyq3Cc56dNBWOp6B+frspbmlVRkBQp
93cR8mvvTOo1h/nor+dcTVa2E3Bn/f37AqvRrx4HhZtjv3EsSc07vIoDPys7
lNdYKYL8O9fPX7BfLDW9yCSAgFcLBSS8t9BAyZzevK6EAR1+EyYOnTQ0Xxav
NPhHQHfWzvMv0f/dZfmFDKuFYdNo6QZh9DdJKyW521sJ4E6bblViUkAk9YwJ
DjFBQidRthD9XK2+5Mcy5Le/v/ZUJutiv14p530O+ev+X9KHWcuGESO3T00R
yI/POR+O4v64BxWfvo35Um6dILcA9akQ5O3sjLx9f15Nae4cBcJXBh7m/hli
qRoNL0pTp6D+MY9aRSQJ2lt/xgYRHHiv5MY2nEeAQAWHGf9//pwsl2ubkwPr
o0825WAfi2ipvHUc+abF0LYiGfvSTOmf0e496BfbA9Mf/aVhbb4HbzP2dXu2
jc69U2zY43y0gW8lAUc0Ol0T0b8Lbp5P6m8hITCvXKIX+99TVuXsnWZZuJ63
O+2tLgXxVX9OX0U/MBtzkigX4sDu0vsTHz/8Zn3fXbp3rRsN5guyc6s6FED1
geyhumgKFi7KjAtAf9SZNN9QFkWDXc9mS4FQ1AczbvbId8zT7Agh+jwJb2zN
5yVNol9ld8T/QJ7YmPzt8BvUh03XDfoSzi9XKVskAr/nbc/+6PXzD8tsY9nK
Ll02XEz0VV9TxoSxCa7Plgtp8NFZRegir8u7pXOJK1Kw6KDWFFsU+9Hm1n+C
fPGwxb/oYImcGPjPXT07cJaGtHfl7CNHCVgR0h6ZjDypXh23fNdzKfi4wC3a
SJoDn/Wq5KrGCdD/FZhjIBgPO2/ZNzVIEjDgR3p0XKNgqO7wgyl8nwvqbj+e
ERSIHftUKxNMwtOUd5EmSWxo26VucGAz8lqa2xemHBtC0uWSptEvZ811G01F
ke9bRS9NYD4eG+l77TPBhtQopW0M7HO7mo/vuXKDBp1Bs4uWPiTcsfz0nUcH
+XzX2We7z39liRoLZs3Dfa/e/+3GQ9RL7OZgm2WNFHR4cFz10Z81nJ/uSudD
PxwYlbQ/wYQ+gYcOKqin2npPNY1HBESfKbHvQf4WGcp0/eNEgHeifsTN7TQY
Znhxak8x4dzbwFOKs2ywzCoZrgiXhcBFhmPeohzYGwVFFw8TEPv75rJtmIdm
qn78Z/G8KwpnljIqOPCBO1dUdj0XrPr82OJIMeZNxMXo/39f5Hn2cAN3KA1A
j9rbW5GQy/rpun8tBziPrAcfNykCVD3ffgvzrCo3oSDbXAGitymEu3yiIJIt
fLePjwmPznzoc9qB+/Rmt6ZVrTQ0V5gYW95BvrXoepQyQsD52vb9Xc8p+CIb
kRndIwtHpsp8w0RoWHM/cdW7MgLWTh9ZuGaCAq6f/WXzkSevhAkZPhGlQJVH
+P2RkyRMbpd5ugP9VHK9fwLXZyaknHfICg3jQPSiefOeYn7MFWyX1UnggOfr
UO2769FvCn1NN6RRENg+7h+5G9/vW+gZ1ZM0/O7oSG7+SYCsxbZpR2HUu2R1
2zHkmed2O0/yrKPAtvZBTjL6lUSHsMSwPgURPTPVsjx4verXKu/we4GkoKRd
M99ZiQanJSzb2OAb3ifCpghgn39sp4XnkWeuvb/TVB664t/6qSEvZd+9vdYZ
+Zr3POiLXEF+Eq+nEpXbWWzBr6t2LMX+N1N06G6EPKj8DefWQJ582h4qdruO
AXfG65v6T2G+71Jx+Dn2hbXkgztDCN/npGJs4x0WCcPfa24J4z5st7SzPFhK
gHuB4+7+VAou/BYy80vHPvJScIEh+j2XlfjF88jHYwU3ExrRP54PNNd4dSrA
OdZV79Q6CrTXdd/y88frSb/4q8ZFQ2PFm+9d2Mc0e1I16rFfdtvpjzXvJ6BN
NdzpIvL1wod5meefSkGoX8DKeGP06/6MeWwBJjz9Zm80dImCB7GatlJ8MiCh
LrD4Nc6ffccrfuMCErLM5/unYR97fGihRAHy4YsL0QPK2DcdT4mtPoi8FbpB
/8vwXTZQenGk7XwClD12ySTy0KCWxnffCM9vQ81Bo5v5NCypbardjjxSem+C
7VjOATcFJedUfN72D17rbcNoOLDRV9Y1mYDZ2DKHMwIcKP5t4JUrhXpPk69J
v0aDZq2I3VtnEtZNZjP///3OgWXVdUcteWHTgYiw73exb7zTmc9YhjwSdbnD
7z0F518skBE1lwJuxosl0w9puMDzwPQV5omnt+EGXZyX7pb3wTyi8nDMy8Ha
FnntQe1w04sMAkaYY1Z7f9BgH/9sTTGB+3WKmOJBXjtVnKP2KpQBOqNJj6zQ
r3csOP3PKYKAkxc4M/kVbNAeBrnaLSSs96uM/YrzaqQfPozH9zdqce74eouC
n/MEOpYRTNTpwnuHefH9FJTWliLP3DZNpX3SsI+++ji9HvUw8zKi8dYwGz6d
E54LlED/M9m3oBI4cFlkbP+GvwpQr/TGtgqv96/3h6vGKwI0V63btjSdDRJW
23NfepBg8vXYcy/cr99nLqR8WkuAc27ck/QYCkzOrR++ukgeYtgzX2X6KCgK
mLoqxcuA6ofzswbXY9/vdboejvkwJaGSMR+/t8rrmTd1iAT/oz+mU0gOjLPX
TO7Fz2adlmJLROPhfd/xXj1HRRCd46p/j32+ToOZVPZEHio4WgL/FDnwJEDQ
7Avq786vHrm6YBoqchpuWOuQ0NhXF1aVzQGRzr4DJ4wJyFe92yN8gIbWn2F2
GwKQV4N5wp3Rj5pTTgiZzTBBePm76tDbNBwO/bi+pgv7v0Sy8Az6z6nHD7JW
ZGtD3EhusGEyBUFzd1fYRhIQ/3XDs5o87Gs6x8I0LEk45JUhcMEG+xK3Z0wk
3l9gy9cYY7w+n7jz/WLkzak0cdpdjwYrfvnfJtgfM4Jez4m+pYDx2q1iSywT
VE2zw9SCaDj9+oaWQps0qAl+ygt/QwN5buVeBu5jo4RthjjmhcjQwohDrxUh
ZVAxXf0PDfHsKWHng70sjXmbP/55Q4HpP1dXMewrXyra6nIPUxB7IYQrGvsI
Kb589DX6Q9HtQ+26cwog76Dj4op9+6rN8m0pWxUhRPXqtzX6HJjecV49Dp/X
Xjlp/10pCvbZP9h747Ii7BJ7tPUo8uSuHU/+bD+K89ApSdd/wIa7W2v9vZVI
GIlckjfjT0HvP86biSJ5WLZkhbop+um2G1t5WlejHlYnWVNDmG/yhvph3T0s
zRTt+81ZNPirskxssL8EnAnOECqhgGPirx1dIwstL2iV3SspuLV4e8wjARK6
9kuMqrciD6gqPDLEvC0OnIJ9HRQEz9uQGxPGBLciLsOLERSQ7oa9Tvh9Y962
Aku8H5XF4+02LQcuyoINytIUNPKcBuOphTBQnTcaHUdB5cULe3etYUIpb+Rx
Kewne2JU9poqyUAjl4pIAvr/hd7YXCPsa63pO3J6t1Aw79aLFatqsb+qNbQs
OULDyh2TPkvkxUCgsvZz3R4OtLbNZEkJkbAt1jtDbwkbPuduZueXKMI9SYH4
oQwapBhhzQLaSvCJn+YEI383B0yIVCO/5RexhJwxb6feG/QNvecHqcJb1yg2
G+7sL790F/V6+ds/P42VNPbh9pj8rA8sTT9SQWiQAtmy73qHA7hgKYMAt33Y
H6o40v2Yz08P3P5ZyaLh08dvTZX//72Lzmz15o0cYBGzgiGbGfByNFOv9wsb
Erm47H6jP3YYzjQ117ChSOtwPyUgBNbpdQ+XYd+NnBFQZbYpwmOV+xVyqTRc
YvRYleuS8F2h/IScMAf0N67z+n0D+2+kuXvBZRretSxuf+fCgKtjh3sPjlPw
TL2z3PmAAmyd+JBZphAPptoM8T70W4ve3iQzTw4oZ9QN3DJA3tmn4HkygAOK
9P4OR8zX7B1aHN17FMD3DZzZ5Zj3Ugm2FhfZIPJddJPNc0lQ75LTtyhgQ0LJ
Ae3wMQLEdFUqlHB+rF91yRln5OGWbwxf00X0X97z93dcVYDTl42VzDH/ilOr
Lx0WHWO16DoUR4tz4N931kkvawKGNn1vblrNwee8tM/vEgHVxR+XVGO+OVVJ
PnbiXQRWXsve1lyiIdTtUcKsGQmPDrZlf8X5hxz+cD7ozS/Wqin2vi+FyF8F
tuXFGuiX30y4ZuQwLxqL1St3kLAyUmB2+28KtmhbuPow5SB8plYqyJqCM7ah
e6RtSHiV9zKzrgP93VPhRR3u7yPFrT9uL8d8j1S6dhp5svWuaLb6AxpGlfhm
7S5NshIDNrbwhbBBf8H1u67YR3Myt3cuQ/8p2LLw3pZtJJyJjDmrhHlhFDC8
IhP7g4rN6okUzKfO1VfiTCT+sRpE3qnb/2OD/YmHnqrYB161+iQ7yHGgyi5m
Xosd5pOUY/q2yxzotXl5uMCEgIwJo2TBRAo5YS7aB/2EsXS6qQD7MX1z7wSb
nwT2QvPHO7EPMa0tYoxjSZjuVgq4gn4zM5z+Og35InunR1aSChuCdprWGT8h
YEz5R8dTPM/zOqrvy48qgPV7ze0hFA2JS3U8Dh4mYaLrQY80C/02M2DrrhAC
Ru9M9oahv1VPW4ccyhpnlV1L9T12jgYThfXtXBZMaJ/3rSjtPvLfylvjo+VM
eMKaWaqL+SX7okFAF/t27D07Zxkf7F/1bWylEgKsD6hZMbH/8Rut6Q75Lg+f
Z0MjvfagXi4V2RMu2P8SA2UsYml4OOCQZ4082OewMafHHb/XuvtMV42E/BM+
izt9aFjVqeIphP61q7hnMnoN9tU0drkHnzQY6YcGiLykYGKJYxQ7UQDmz+32
eqTGgXMFJi47MI8/vjrjWY38pfl5GfkX85sZMfL0eQoFUhG6Tb/tSaBdD7rn
K1Hw6H35Yjt5EiJe8i+NysXzSVjdexL57jljVfPHTArm/6zzv3CNCR99naTM
LyCv6a1T9cc+k+s767II70c8fzoS1UqAcOijQ5PI18K2f6N2Yh+XtMsIyEW9
GVw6pX8L92WV8smdAok0+D7fatx5i4Cd0aV3aF4KtGa+1zccI6Fzk5Sh7g0K
rs2voFbXMOHYWi+FXuSHYzm9zNvII7eXy0yketCQtyRhrvoLAbZHawv3vKXB
2P20VLkGE+we8WkSUhx4czR9oY6sJCzx9R/WPMABvprKxxO430ZDLC8f3E+f
spS6VcjbZy595Q7APrfuoT57/NsMayp4gliRwoZ/faE8MlLy4JftxuN2D/1X
f+eDk2NyMOfYydPRyYZfMzvuXf9IQPjdDa+2V9JQwuEx9gokgDwT2Nb3gwKZ
1m9M3RcyMHCdPS7rRsHWNzVSp/SkYJNy0mK1fRwIPpG2+kc7Ad1n5m9h9aK+
vn5s/+DKhCaP9jOvsN9MS3tpu+D8vderRlXiPq3fEeJ767A0PPHsefsyEPVq
IO50wIoX4u/nHBK4SYPDYxNd0aMM4KlekDj8gAOHj9JGeZW4T9XSktLIO5Fe
9fzqrxmwymLr8vceHKiN8XHoMGPClsArz4b88H4H5SzGvUUgewX/Bg9vfN+7
LF3upwzQ8FmUysbzWHLZ1rgb9fLhTG3MChMaZs7Lv+BCXnbV6DbklNAwoqpW
fM91kBWmNbFry/+/z9B7+etPNPL7qX7rUgEKXv699/dskjS0Zzx6rz5DgVEH
z2HxpSQ0jZ2adDtNgf+Pz+By6icr5NG0Pfc3NpSI8D7ysmRCMa12NGKAAuN5
j07v0SKghb7KaDKkQdDXLCrJjQTlD+BmNsWG3LUvjq5eMc7q5ETJZ/jh/DJU
u4LeysN+B53xAOR1M6m6H1LoJ2cMv/Rd4ca+0sbZtV9JCkbIlc/lCil4vX/f
IrUjBKQ27AjvQN6VKF3v9SpECJQ95V+tycD8d34frKzDAJrfs+fCEvQDg01X
xtKZYML3Z81t1KMOl0b+c9TzpbJrF73/USC3MFPtgIEYyP84uF0F+e2a3GBE
wQkGiOVWhvuh/uOHCWH+kAGWduHt/H5ZCkZ23nu8ZvUY68Mb/h9hG9nwZCOj
ySCchDZv9fbD+2nY1C1ygV9QEvY08/ZsxXzlqrIf8Xb5yWoouKexu4YCR9HF
y9ZKk3Dkzr6B98h7b9daS7lnKMBDmYtHTzcgzxpXdx01ZEDJE1u3pikKPKqU
xOJlf7MUA9l2P4Y5IOD7erzVHfOvUOc7h0FBkxd4HEvE99mlzL2pCXnIVXj7
kcGFYHC+dH+HLz6/6Qdnw6NdrHThe1m+dWx4JZTkMII8t24wi5LU4oC5wLjF
V3NF0J47bqt4gYbLzZyot50McHN6G7T2Ce6f/AtiThj93pQxpt1Kw4/sUo/Z
MEkosPieNIt6fVd0X3rrM0VYS+aCxFMaVK7DM+tRKVBccNr+eCYNeoYqz8/m
8UPB+iy3vchvP7U/BJ44qwBMvWbJxcgjQ9fkzig28kOZyFZWNr7vpS9Pty0T
54Uaq11VMbh/3PlLRtQxD03ZXxdfyGLDuqNPHUKdkJ8v/Y551s8G7n0z1oA8
nNXpVDeA/kTNGH0aqycgUcq/YvUwBemrdOp8cP+2NHzgGWWiX77OOV60igS1
RflCPNYcOHHr5EWjEC7w7KVLV5dR4PK+yf4F8pJQsahl4hs2jAe/FvZ2xX5c
eNI2uIUD2uYHG+fQj86HDuzUNse8YG0frl/CAIl8c5vaPxRs3n1yn6kjCdW3
w+8VpdDg1ae/xNOXhIalSVN/XThwRcYtyC5uktUsnf2srY2Gm37P5+rEmUAm
eremP0aemogRJ5HXTW9Xy1p9jYOcxBDZhT4yoFG3kWhBv/z+plUtCHlzdYtl
4DT6nVOu0ZTPDwJz/bPe0QIKVD+zrqbwkhASUJd7WZuC6YpUdwnkk5fbEuQK
kJds1dO3+OL+2R/KuGmHfffc6Av/VWO/WAVqh4y6kZ//iv10tLzHhDPBlUu3
36VAcImu+VXUf1HGln9DZ5F/Bw8U8uL1lTkS36bFaJDnf2qrv5gE3lkv4NGi
4bXrIROZR/KQKrZbXxj7ir3VYpXHgz9YoeK/QwnsP8cdvpVNOC2C1Oaq6Xg3
GuYMtNyMkfdWOalu3Y75uOVyt6PbKPaRm3f22OD9qvx2N8xVMUHx/vC1jQ85
EOhCiUojP6y9YORxq4aGLo7o6HHs3w3/Poe0q2A/UtHm8Y2YZO1xfRaRPo8G
53uv9UwxH/39R6OuhGKferf0kLAeE55ynnA2FOD+7Qs6FyyOvCt95+FD1O+S
naufHvkkBPlm3yo0v9Aw4VMq04Z5puW82bpehoa3V8946HpKg+BLM8uDOE9b
L/5N3vj+BXRkaowzBTtzT9gvnJCFmRhHxweBFPy68tH7UDETdKHQstiChiJ9
vc4LyE+5ubZfKVMK0kadYnqDmLCb2qPUh345Qscdk0be9XqW4CswSsFNgx+B
bm5MUPa5qijazobNG8UOfnCdZNl/1ROPRR5u5jdKB0kFCFYXPLrFiA2Pub9M
LNMUB1rf+95JA+yjryTPVv3P81VL1t5WwP6VHT2h2zrGuutixbUvDvlxmcbX
UezXb2sGmk4hn/5ZIR8szsJ97svmEe6mYPWkh8KbQQaY6RdntDRTkMuyX5i5
hIRwtXMD18woKBwjZJTUpSF+rNQsGPM/a9e8ZnvsJ8U27+Qe4j4PcZrKDfjw
s+A1/pAcNjAfrlbgSJKQF3qIxdNFQdbGBi6xggmWyLUoD9NdFIQlcF/b/5WA
F5tnFxx/SYPwz0+tpyrnWGfoQ5df7sL5nBI/LrGdhMCxurw8Hho8D50IuoI8
fSgjICkSeambX0grspCA4ti0pst2FNgwzt19pcmECXFobdqNedJxouqoDwHH
7ycFcrDvvo5dovzzFQGe+rv/+ShxgJ43mt6/i4TDydbqxdgH+xZeE2jE/fEd
lzRpa6Hgh3mxPmOMAYJxDywF9WlQkOsI7L1LwEH/+/Wqo2xIjYx5+i9qgnVV
rswlGPk2LUavYa0uA1pHI/iOiNAg80WyS7BNFnn0SO891G8wT93S2ClpkNjW
+0gBn+fPy0qVSQsCvGY/b1i6igaLVdJC8cjj7Jc5DnLot8/UZib8oghwbsxM
DF/MgQZPurNfEc/zpgJxYxT1YPp6XY4sCTEFUg79+L4n3B8IC3u+ZvU07Ti6
CPnbw2WGl5cgoeXuC69gfg5s9nnklIx9ymLcxkYHn6+728JEcoqA+xEgr1CP
P3+zxn2PNgm16m+ecHmh3l9GSGqI/mA1uAz+7kbe3vvGRmHdegK2D4oxHLCP
GvU5PVDaN84a49c/0fSEhtUGvuvWuzLA5bnRG9KSA50KJzv34/2n2Tn1oqsp
mJXNquwWk4HFHtcduDD/HevfLK78/++bBKvOkP4cWDk8GUE/UoTYB958JZg/
N9yiKloFF4Gl6d4b785SwObl5YpHHpW94ne4tRPnaS/8Ot6WgG0v3zk++oR5
6mR77pIVAdnVyWkeoRREO7+//ODjO1a6ZsP+85j/zUeHNdnon+POs2GXt9Bw
0GnGsOOAPCTdNFwig/s+z9ysfO8GAsyuZzKC+dCvfxaKXcvB681vWRVgScPi
9sNOTOybyhsW8NegX71kC/tmol/5xw6pfv1Fw8c19/2/8BDgvu/fijdEPPCb
3PppaiQPgsniYUpcNATJW9AdBBNc/vL0G+phfi5zaUzG/inrYKzPeY98cdW/
XEsCeY6sgL3lNCiddlJ6g3yoSbweIpBnIp9+c79wngCZu/5Riaso8O3Vy/Ff
hH3G2PZHG77vLvGRE/X1suAa1qb44hsNR9WNy4bqGfCI7+U7F5y/3lY34wId
EoQvics1tlJg9+RKh54yCS5PpI+MZLNhQVTI41j8/1rHq+LzdlKg1m1X/A77
44e3YdfaqtCf3PfGMoTFQOzF0472rxSEH80sYKJ+prfd+MhRjIcjcuFeM+gX
7n2vpXxaaLjC79q/biX2z8pVJm8ICh5L+yyyOIx9vSoufc4LeWngsvrNhSQ8
sDUKcEE/TPsqnSSsiPtloLd6FvMoJv+kv58MAQU3qeMTwTQkhZ85ewD9ZJFe
tLsv7lOfu5yMJ/KAWs3mKcdsCl5tFNaODyAhYdeHgKEXyOPz1K+E5LSxuP3z
A0KvcyA5LteG25OA/OaShoZXbLguNCHq/fM3S+eOr2jm/38fxHKhPK8woeWW
kl6AFRsCe7itl6IeTO047D24D/q7Y4RsdhKw22Hw2zka/XRc4JfGERIq2z7n
aA7RsKDj+8A/XyY4WzX7fPCkIYTvbaEXNwkOTj1fvxZRYLLXY+FEjhj0t1k7
FzI48KBDTdVITx5in2eGaDEpUEnJ073uIwaJK7ml1xynwb1ac8Ia9ZGvFr73
/Dbk1dMmxo5N6EdanYtH4ilw9jaRoNAPLjk6yTwTo6Cg/J5XBYnzWiUwVTRG
QYbGMr5HHiTs5n+XzcDPY+7zowbzP7Fenhjv+r6ZA00lQ2k6uA/Davt37vOn
AX4Vtgqrk7C8kjO08xgHjLXZ+yKi0e+/SdCGmIe8d27z52NezZU/ZAmjXxGv
lbmtV8nCrSTvrsEzOG+T1kBn3N926b8PA8/T0NEV6HDZQwnOGmT0SGL+KYgK
e8lkMmE6vLymxpIN5sW56fNeE+CnMGl2BXnbIHdqRw72gbyo3BwLEwqmBq3C
S0+Ns/b5q7rnGSCv2R+u8sM+zPs9+wl3BA39bAnwDSHgt02Diu58Gmrzq7as
uaQIM8obk27yUFAiochrPk7AU8O6KGILBT5bts8tR74LXFRg4tuJ1698pesZ
yQ9rjylJ3KmjoMErQmjVWAeLPgBXle5yQOacxKL1Ddgvio9dC0L/S/kcss9X
CXmT8Oi0l+eA9+QCgZZxEuotv3CMPSjQ25DtkYD7aEDS8VY2HBi2lsrqVyFB
ILCSuTmahluFv1xu+KEe9DYwTuI+HV4dKxDXS0DFtIAbbzYH5p3RTyiOQz7M
exYZh31T+vde7qLFsvBqLMx7BvntbLpxqIYZAU0qIfuDPrEhP357jS32Z9Wm
rIi1X9mQsJTwFEDemZLqPhV2mALwPc6KRh41+T63d806DoRt2bG4ZxL3Pzo5
XAh5jK/zmkssPk9SiMzrdw14vufrm9fZkfDqLadx6XwOKDacvuGrqwIW6xxb
VPH8zy2KbC1jkhC9hbMvFvuf0OWRAr0VBFzgSe5V1qZBdxhObrpHgM3qsLD/
rydsdUqjGK+/+0uwqh3yotL3P/SqGwy41RJvdLknDmyeK/7rfsqEeKnwvQqY
x6nKL7M3iJCg5DB6pwj5QqTieXhJEQE3aL7LpsoU5Ic7qRz8/+9rlA/JbRan
YeGuFxI+5gSUK2o6fsb+rnvfunPdKVn45Oyh48NLw4s1UjvTbhDQvJwnyPkI
G/zGirxMaWmouTm8nm1Owyu74657Gbiva3rnzVtHgeW/KLlhNUW4OfJr98Bn
7NsJCWsFBglQSV3DXb4gHq7Ebt6ZGzHHEg29OPBOOB52Dpb5P13DAD4jY+9m
ZRpKOw826GmRULrcqWLbCA33uZm+1Nwgq0UyPuhuCQV3dghsaktpZ51bX1Uk
voiG8i9CeWOLSKBTppe/wf1UGpOJqZaQhqUppxorJLAflUe33sT+MUhWjpGv
sd9JqjqeQH7Lersuv6mIA6cNrdYsfceEGrUR5/FoCub+5fzaqUFCyrWy5XPI
65/WpC7X+t+/Mtlw+jcbTgjueXsQ96Fc6O/przhfuYINkYNmDOid3tu8IZ8N
F7jsbAwP434eEtU61s0Gr65NJ74IkSA6z107Hf0pe94OgThR7Es7DjdtXcmB
nyXrV2leloWTimU572rZ0B6uciN5BQnpx/NaqjyQN7pmRlVMSZhMtev7G0KB
jkHcQfdgEsizk/btqDd2v5mvBep55cD9ogLkIRAo1FHYRoLmJe7AXmkaSm5e
bl46KAPXNR9wGBo0DMY61qvPiUPHzvc2B2vYUOGZyvtcWwRa0ztfS7vg+Yjv
DX2Gfc8ypaFf8wXmjfnVD9ciCPBRea0fK4X5EeiXfBP52Gr76pVxTqhHrYxh
EudvN94cKXWKDR09Omf7nv1lbVup+PNZIuYnJVO2MFkc87Eu5V0YDeEk09AU
+6nGK7s6sSAamuHV6K8rcnD9kTN5DvNnavF+iVrknbPj/L8W4PW5b347eA/7
VOhL1/5H85CntVS076Gflz6fLsuoi4OQsbJ3vqeRZ//pPWeNUCDvPGu7xRL1
ssbSOyoQ+xLz00vVbBKYixnKatdoSNhocq4Z+8SSD0tvJQ5SkPSiVFd+kxIU
XEt3kv1Lg8/k74DVf/lBO+uk5/P7NGQFqLpE4nnk/uk3sViC/BXzNZX3sQzU
p0tkdrE5cEn/8vSu+4MsVydXWQVR9Jsw65enBplw5ZJzWMF5NiSJnncOwPMa
uHlq27YcGnLkHYyjgxShO2Zo0+hb/HyQtb7oEwGrpuxKDFUwj+wqfuRi/m4i
D+h5Ih/IMubbn9hLws9yy7QMK/TH+kX7ktBfvL1s+HmR5741n/VX1paD/LTy
1jnkTbsfntoX8Pxy229uyNlKQ+AVmW6TayRsS4IVfA84oKkj/bXyPAlGlneH
ZrCP56onKyzGPqbHq7SyPZkGt4jUCXPMW9HvbaqNJchrmxWFtLFvSYUE2rM3
Yh6MLl/HmCEgPnu/qf157A8fPb4+dWHAiaVjxV9UabDOcfGRwn6hd55/mWEQ
5sHa8ELePh74D+s8uPc=
      "]]}, {}}, 
   {RGBColor[1, 0, 0], LineBox[{{800, 0.27}, {1200, 0.27}}]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->None,
  AxesOrigin->{890., 0},
  Frame->True,
  FrameLabel->{
    FormBox["\"Steady RNA\"", TraditionalForm], 
    FormBox["\"Max deviation\"", TraditionalForm]},
  PlotRange->{All, All},
  PlotRangeClipping->True,
  PlotRangePadding->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.5168777067519345`*^9, 3.516877733260451*^9}, {
   3.5168778510491877`*^9, 3.5168778648509774`*^9}, 3.516950381953125*^9, 
   3.5171965356875*^9, 3.524383713046875*^9, 3.524398091359375*^9, 
   3.52490209325*^9}]
}, Open  ]]
},
WindowSize->{1909, 952},
WindowMargins->{{Automatic, 0}, {Automatic, 0}},
Magnification:>FEPrivate`If[
  FEPrivate`Equal[FEPrivate`$VersionNumber, 6.], 1.25, 1.25 Inherited],
FrontEndVersion->"8.0 for Microsoft Windows (64-bit) (February 23, 2011)",
StyleDefinitions->FrontEnd`FileName[{"Creative"}, "NaturalColor.nb", 
  CharacterEncoding -> "WindowsCyrillic"]
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[545, 20, 485, 9, 70, "Text"],
Cell[1033, 31, 135, 1, 43, "Text"],
Cell[1171, 34, 7031, 179, 538, "Input"],
Cell[8205, 215, 238, 5, 43, "Text"],
Cell[8446, 222, 6143, 135, 794, "Input"],
Cell[14592, 359, 190, 2, 43, "Text"],
Cell[14785, 363, 1221, 33, 174, "Input"],
Cell[16009, 398, 333, 6, 43, "Text"],
Cell[16345, 406, 7564, 179, 491, "Input"],
Cell[CellGroupData[{
Cell[23934, 589, 135, 1, 90, "Input"],
Cell[24072, 592, 108, 1, 60, "Output"],
Cell[24183, 595, 3567, 63, 250, "Output"]
}, Open  ]],
Cell[27765, 661, 98, 1, 43, "Text"],
Cell[27866, 664, 3267, 77, 258, "Input"],
Cell[31136, 743, 542, 8, 43, "Text"],
Cell[CellGroupData[{
Cell[31703, 755, 12694, 297, 1048, "Input"],
Cell[44400, 1054, 91, 1, 33, "Print"]
}, Open  ]],
Cell[44506, 1058, 104, 1, 43, "Text"],
Cell[CellGroupData[{
Cell[44635, 1063, 97, 1, 61, "Input"],
Cell[44735, 1066, 137, 3, 60, "Output"]
}, Open  ]],
Cell[44887, 1072, 88, 1, 43, "Text"],
Cell[44978, 1075, 2212, 47, 230, "Input"],
Cell[47193, 1124, 85, 1, 43, "Text"],
Cell[CellGroupData[{
Cell[47303, 1129, 199, 5, 61, "Input"],
Cell[47505, 1136, 204, 3, 60, "Output"]
}, Open  ]],
Cell[47724, 1142, 105, 1, 43, "Text"],
Cell[CellGroupData[{
Cell[47854, 1147, 967, 23, 61, "Input"],
Cell[48824, 1172, 189, 3, 60, "Output"]
}, Open  ]],
Cell[49028, 1178, 117, 1, 43, "Text"],
Cell[CellGroupData[{
Cell[49170, 1183, 529, 14, 61, "Input"],
Cell[49702, 1199, 11344, 344, 403, "Output"]
}, Open  ]],
Cell[61061, 1546, 232, 3, 43, "Text"],
Cell[CellGroupData[{
Cell[61318, 1553, 738, 16, 61, "Input"],
Cell[62059, 1571, 210, 3, 60, "Output"]
}, Open  ]],
Cell[62284, 1577, 225, 5, 43, "Text"],
Cell[CellGroupData[{
Cell[62534, 1586, 890, 21, 61, "Input"],
Cell[63427, 1609, 16593, 506, 403, "Output"]
}, Open  ]],
Cell[80035, 2118, 225, 3, 43, "Text"],
Cell[CellGroupData[{
Cell[80285, 2125, 844, 18, 61, "Input"],
Cell[81132, 2145, 230, 3, 60, "Output"]
}, Open  ]],
Cell[81377, 2151, 115, 1, 43, "Text"],
Cell[CellGroupData[{
Cell[81517, 2156, 1006, 24, 61, "Input"],
Cell[82526, 2182, 14560, 445, 403, "Output"]
}, Open  ]],
Cell[97101, 2630, 164, 3, 43, "Text"],
Cell[CellGroupData[{
Cell[97290, 2637, 1598, 44, 118, "Input"],
Cell[98891, 2683, 20954, 350, 405, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
